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Abstract 

Entropy is generated in high-multiplicity events by a dynamical sep- 
aration of strongly interacting systems into partons and unobservable 
environment modes (almost constant field configurations) due to con- 
finement. The effect is demonstrated in a non-relativistic single-particle 
model and a scalar field theory, where it amounts to quantum field 
Brownian motion. We analyze the quantum decoherence of partons, 
which formally corresponds to non-unitary time-evolution and causes 
entropy production, in terms of Schmidt and pointer states in the non- 
relativistic case. For the coupled scalar fields (partons and environ- 
ment) we derive the Cornwall- Jackiw-Tomboulis effective action and 
equations of motion with the non-perturbative time-dependent varia- 
tional method in TDHF approximation (to study the time-evolution 
of model structure functions in the sequel). We obtain a model- 
independent lower bound for the entropy in terms of two-point Wight- 
man functions. 
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1. Introduction 

The idea of highly excited hadronic matter can be traced back to the earliest 
observations of high-energy cosmic ray events with an associated high multiplicity 
of secondary produced hadrons. Most notably Fermi and Landau in the early 1950's 
initiated a statistical approach coupled with a hydrodynamical description of the high 
energy density matter, which is formed and subsequently evolves in such prototype 
strong interactions [[I], 0] (see also || for a personal historical account and further 
references). Until the present day their pioneering work has influenced further devel- 
opments in this field, even if for about two decades there existed no hint whatsoever 
to the most relevant if not fundamental degrees of freedom of high energy density 
matter. 

The statistical approach, in particular, is motivated by considering the large 
amount of initial kinetic energy which is carried into the reaction by the colliding 
particles and is ultimately transformed into the large multiplicity of secondary par- 
ticles, i.e. produced matter which finally flies apart. By a fast (on the scale of fm/c) 
compression of a relativistically large amount of energy (^ lGeV/nucleon) into a 
small volume (on the scale of fm 3 ) one couples very effectively to strongly interacting 
degrees of freedom, i.e. those which materialize as the observed hadronic particles. 
Loosely speaking, a large number of initially hidden underlying field degrees of free- 
dom must have been excited during this transformation of the apparently simple 
(ordered) initial state into the obviously rather complicated (disordered) final state 
of the collision: It seems as if an enormous amount of "entropy" can be "created" here, 
which experimentally increases with energy and complexity of the colliding hadrons 
or nuclei. Where does this entropy come from which is one of the characteristic fea- 
tures of multiparticle production events? How can it be interpreted microscopically 
in terms of a field theory of strong interactions? 

Further on, working towards a statistical theory of high energy density matter, 
Hagedorn led the important intermediate step of taking into account the informa- 
tion about the dynamics available in the 1960's in the statistical bootstrap model in 
the form of an exponential hadronic mass spectrum [[|, ||. Here, as with all former 
and later models and their phenomenologically rather successful applications, the sta- 
tistical approach immediately suggests the crucially simplifying assumption of local 
thermal (kinetic and possibly chemical) equilibrium among the relevant constituents 
of matter in very inelastic high-energy reactions. After the advent of QCD the cor- 
rect theory of strong interactions is nowadays believed to be firmly established with 
quark and gluon fields as the fundamental degrees of freedom. Thus, corroborated 
by the observation that asymptotically at short distances the running coupling of 
QCD vanishes and once again guided by the equilibrium assumption, it has been an 
ongoing and ever stronger effort to identify theoretically unambiguous characteristics 
and find experimentally signals of a hot and dense "quark-gluon plasma" of almost 
non-interacting quarks and gluons. This intermediate new state of matter should be 
formed at sufficiently high energy density during hadronic or nuclear collisions (see 
references Q and |J, respectively, for early and latest overviews of this field). Until 
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very recently, however, the equilibrium assumption ( "thermalization" ) has been com- 
pletely ad hoc and is still lacking any detailed justification or deeper understanding. 
Even if many global features of high-energy reactions can be reproduced in thermal 
models of varying degree of sophistication, it has remained quite a mystery, why they 
work so well. 

Connected to the problem of large entropy production mentioned above, the ques- 
tions arising in the context of thermalization are: Why does the thermalization, which 
implies a saturation of entropy production, proceed so fast (on the scale of fm/c)? 
Which are the relevant degrees of freedom in terms of QCD and with what kind of 
effective interactions ("hard" perturbative vs. "soft" non-perturbative physics) that 
contribute most efficiently to thermalization? How do these "randomizing" degrees 
of freedom hide again during the final hadronization process? 

Questions pertaining directly to the latter hadronization or confinement problem 
presumably will have to stay off-stage for a while. However, even those concerning 
(time scales of) entropy production and thermalization cannot yet be seriously ad- 
dressed. As we explain in this paper, the necessary quantitative studies have not been 
clearly defined from a conceptual point of view. It is our purpose here to present an 
analysis of the entropy production in very inelastic high-energy reactions with an eye 
towards the parton picture, to provide a plausible hypothesis about the underlying 
randomization process, and to outline some calculations under simplifying assump- 
tions, which should strengthen the proposed scheme. 

The plan of our paper is as follows. In Sec. 2 we introduce the definitions and 
formal apparatus which will be applied to analyze the entropy problem. We employ 
density matrices and, especially, we explain the notions of Schmidt and pointer states 
in the Hilbert space of a complex system and of environment-induced quantum de- 
coherence. They play an important role in studies of the measurement process and 
the classical limit of quantum theory (Ref. j?| provides a useful introduction with 
numerous references to original work; see also the review ||). Recent attempts to 
formulate a consistent approach to quantum cosmology initiated by Gell-Mann and 
Hartle and others have given a new impetus to study these questions again [[7], [|, 
which touch upon the foundations of quantum theory. Basically, there one tries to 
understand, why the universe considered as a closed quantum system (in an overall 
pure state) looks as classical as it does. 

For example, low-energy propagating modes could be coupled to dislocalized dis- 
crete modes left over from a topological phase in string theory, which latter ones are 
hidden from a low-energy observer performing "measurements" by local scattering 
methods [[L0|. Thus, there is a dynamical separation of modes according to whether 
they are observable in the low-energy limit of string theory or not. The unobservable 
modes naturally constitute an environment which may induce quantum decoherence 
in the subsystem of the observable modes. This in turn leads to non- equilibrium 
evolution, i.e. non-unitary quantum mechanics in the subsystem with an associated 
direction of time and entropy production, in particular. 

Analogously we may ask, why a pure-state high-energy scattering experiment 
lends itself to a statistical description characterized by a large apparent entropy and 
possibly quasi-classical propagation of partons. In this more down-to-earth context of 



2 



strong interactions similar phenomena should occur on a microscopical scale (Aqcd ~ 
200MeV) which can be related to a dynamical separation of a complex system into an 
observable subsystem and an unobservable environment. Surprisingly, the subject of 
our present work has not been analyzed from this appealing point of view up to now. 
The paradigm of confinement has been around for about twenty years and it seems 
unavoidable by now to think of unobservable quark and gluon states or, rather, long- 
wavelength field modes feeling or producing confinement as constituting a dynamically 
active environment leading to quantum decoherence of the "almost observable" partons 
(in the sense of parton-hadron duality or deep-inelastic scattering), to which we have 
become used. Thus, the situation is opposite to the one in string theory; in QCD 
almost constant low-energy field modes constitute the environment for the "hard" 
partons. 

This may have important consequences for the parton picture as applied to nu- 
clear collisions or complex reactions, in general, where factorization and assumptions 
about the quasi-classical propagation of partons between successive hard scatterings 
come into play O, T2\. The so-called "sudden approximation" underlying the parton 
model, which can consistently be justified within perturbation theory, could fail or, 
rather, be justified on a deeper level precisely if and when the strongly coupled envi- 
ronment becomes important. After all, the commonly referred to analogy between a 
parton and a hard-struck electron in a crystal can be rather misleading, since only in 
the latter case all relevant interactions are weak. 

Recently there have been speculations that so far unexplained decohering effects 
on the initial state multi-parton wave functions of colliding hadrons or nuclei should 
contribute considerably to entropy production in these reactions, see for example 
TBI, KM. We show that the mechanism for entropy production in strong interactions 



at high energy is hidden in the active role of the environment. 

In Sec. 3 we present a non-relativistic toy model of a parton linearly coupled to 
gluonic oscillators which serves to illustrate the formalism of Sec. 2 with environment- 
induced quantum decoherence and entropy production in the partonic subsystem. 
This model bears some resemblance to a non-relativistic electron interacting with 
the quantized electromagnetic field. The oscillator spectral density of the gluonic 
environment, however, is by construction rather different and receives its dominant 
contribution in the infrared, in particular. The encouraging results obtained in a 
short-time strong-coupling approximation employing the Feynman- Vernon influence 
functional technique will guide us in Sec. 4, where we begin to study the problem of 
strong interactions in the parton picture. 

We propose a relativistic field theory generalizing the toy model of Sec. 3 consist- 
ing of two non-linearly but locally coupled scalar fields, which we assume to represent 
the partons and the confined (or confining vacuum) environment modes, respectively. 
The purpose of this model, which may still be far removed from a phenomenologically 
viable representation of QCD, is to tackle the particular technical problems expected 
in the field theory context. Loosely speaking one could also say that we extend here 
the study of quantum Brownian motion into the realm of field theory. As it turns out, 
a functional Schrodinger picture approach based on the time-dependent variational 
principle of Dirac seems very promising WB, [IB). It allows to transform the usual 
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Heisenberg operator approach to applications of quantum field theory, i.e. mainly 
time-independent scattering problems, into a set of manageable variational equations 
describing the relevant initial value problem. Furthermore, an important advantage 
of this procedure here is not to be limited to the study of weakly interacting systems, 
as it is the case with ordinary perturbation theory. 

In a time-dependent Hartree-Fock type approximation we can calculate the en- 
tropy production for any scalar parton model field theory completely in terms of 
two-point Wightman or correlation functions, which seems to be a rather general re- 
sult allowing interesting conclusions concerning the real QCD problem. In order to 
learn about the time-evolution of the two-point functions, of course, one has to specify 
the interactions in detail. Similar in spirit to the BCS theory of superconductivity, 
where it was realized that particularly in the strong-coupling limit the "action" is 
around the Fermi surface in a conductor ||17|| , we propose here to distinguish between 
modes which are close to constant fields (environment modes) and those which are 
not (partons). Thus, a Momentum Space Mode Separation is the essential feature of 
our toy model field theory. 

In a sequel to the present work we further evaluate the model of Sec. 4 and 
explicitly calculate the time-evolution of the parton density (functional) matrix for 
various physical initial conditions. It is important then to derive the relation with 
non-perturbative (spin-, colour-, and flavour-averaged) QCD parton structure func- 



tions [11, [12|, y|, which are experimentally accessible. Furthermore, we will refer 
to the parton cascade approach for hadronic or nuclear collisions ||18|| . This model, 
which is based on perturbative QCD to the largest extent possible at present and 
singled out among existing "Monte Carlo" procedures by most clearly specifying the 
model assumptions about non-perturbative soft processes, is particularly useful as 
a reference frame in trying to understand essential aspects of entropy production in 
high-energy reactions. New effects implied by a dynamical treatment of the unobserv- 
able environment modes, as initiated here, have to be contrasted with such a state 
of the art description, either to better justify or to modify it where necessary. In 
particular, assumptions about the factorization of multiple hard or semi-hard parton 



scatterings, i.e. quantum mechanically independent rescattering |TT], [12], |i8| , deserve 
a closer look. 

Finally, we have to go over from a somewhat ad hoc model to full QCD. We pro- 
pose an explanation of the underlying decoherence/randomization process in terms of 
unobservable soft gluon background fields, neglecting potentially important soft quark 
modes (e.g. chiral condensate) for simplicity. Presumably, this necessitates model as- 
sumptions about the properties of the background fields, since a non-perturbative 
ab initio treatment of QCD seems still out of reach. However, detailed and phe- 
nomenologically satisfactory approximations are available here, e.g. in the form of a 



"stochastic vacuum" model (see e.g. Refs. |T9[ and references therein). They can be 
implemented by suitably modifying the method of background field quantization (see 
e.g. f2(| and references therein). 

The studies and further remarks presented here should indicate the next steps 
to verify our hypothesis about the nature and importance of environment-induced 
quantum decoherence in strong interactions, its particular relation to confinement 
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and unobservable long- wavelength QCD fields, and the ensuing solution of the long- 
standing "entropy puzzle" in high-energy reactions. 

In the final Sec. 5 we discuss our main results. We also point out further inter- 
esting applications of the formalism developed presently. - Throughout this paper we 
work with units such that h — c — kg = 1- 



2. Classical Behaviour and Entropy in Pure-state Quantum Systems 

To begin with, we want to derive the Schmidt decomposition for the density matrix 
of a closed quantum system.^] - We may think of colliding hadrons or nuclei which are 
sufficiently separated from any other strongly interacting matter as an approximate 
realisation of such a system (neglecting any long-range interactions) or the universe as 
a whole, which most fascinating example, however, is presently understood the least. 
- Thus, the complex system is described by a pure quantum state which evolves 
according to some still to be specified dynamics in an embedding Hilbert space 7i. 

The corresponding density matrix is given by a hermitian operator, 

p = l*X*l , (i) 

which is trivially diagonal. Suppose now that we are experimentally limited or, rather, 
we are forced by the dynamics of the complex physical system to observe only certain 
components of the complete state vector I.e., we deliberately consider a fac- 

torization of the Hilbert space into two predetermined subspaces with orthonormal 

bases V = {| )} ("partonic subsystem") and Q = {|})} ("gluonic environment"), 
V 

respectively.^] Then, with TC = V <8> G, 

l*> = J2 c p 9 \p)\9) , (2) 

p,g 

and, 

p = 12 C P9 c * P 'g' \p)\g)(g'\(p'\ > ( 3 ) 

P,9,P f ,9' 

are corresponding general expansions. Furthermore, defining 

\ P g) = Y, C P9 \9) , ( 4 ) 

9 

we obtain 

l*> = Eb>bo> • (5) 

p 

3 The derivation is given here to make our work selfcontained and to introduce necessary definitions 
and the notation; see also Ref. j21| for an earlier presentation in a different context. 

4 We choose these completely arbitrary labels in purpose to prepare the reader for Sees. 3 and 4 
and will be interested particularly in the behaviour of the partonic subsystem later on. 
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The linear combinations \pg) are not necessarily orthogonal or normalized. However, 
let us introduce the density operator for the partonic subsystem, 

pp = Tig p = Yl C pg C * P 'g \p)(p'\ > ( 6 ) 
p,g,p' 

by or t honor mality of the (/-basis. Next, we assume that we chose the orthonormal 
basis of V to be composed of the eigenvectors of pp. Then, by its orthonormality we 
find 

Const -V = (p\pp\p'} = J2 C P9 C * P 'g = C P- S PP' , ( 7 ) 

9 

with a positive eigenvalue c p = J2 g \ c P g\ 2 - Similarly, we obtain from eq. flU) using eq. 

(Pg\Pe) = J2 c P9 c * P 'g = c p- S pp' ■ ( 8 ) 

9 

Thus, having chosen an eigenvector (of pp) basis of V, the must be orthogonal 

as well. Introducing the density operator of the gluonic environment and using eq. 

p g = Tip p = c P9 c * P g> \9)(g'\ = J2\pg)(ps\ > ( 9 ) 

g,p,g' p 

we also see that the \pg)'s are the eigenvectors of pg with eigenvalues (pg\pg) = c p , 
by eq. Thus, surprisingly pg and pp have identical non-zero eigenvalues, the 

number of which is determined by the dimension of the smaller of the two subspaces 
V, Q C Tt or the rank of the matrix {c pq } in eq. (Q). Normalizing the \pg)'s, 

\Pg) = c p 1/2 \ P g) , (10) 
we obtain the Schmidt decomposition P2| of the complete state from eq. (Bf): 

p 

in terms of orthonormal "Schmidt states" {\p)}, which span V, and { |Pg) } Q Q ■ They 
represent the sets of eigenvectors of pp and pg, respectively. Of course, the roles of 
V and Q can formally be exchanged in the construction. Finally, eqs. (MM) yield 



Pv = J2 c p\p)(p\ , (12) 

p 

and eqs. (g[TO|) 

= Y, c p\Ps)(Ps\ ■ (13) 
p 

Equations ( [TT] ) - (|T^) are the essential results of the Schmidt decomposition procedure. 

To appreciate the importance of the Schmidt decomposition, we proceed in several 
steps. First of all, we may generally assume a normalized complex system state \§>). 
This implies 

1 = (tf|tf> = J2 C P = ^VPV = ^gpg , (14) 

P 
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where we used eqs. (11) - (|T^) and orthonormality of the Schmidt states. The sum 



rule ( |14D suggests to interpret the positive expansion coefficients c p in eqs. ( p!2| . [i3| ) 
as probabilities to find e.g. the partonic subsystem in the respective states {\p)} p3 |. 
Such an interpretation is confirmed by calculating, for example, the probability P(p) 
to find the partonic subsystem in a state \p) G V, 

p(p) = (p\p P \p) = Y, c p\(p\p)\ 2 > ( 15 ) 

which is equal to cp if \p) is one of the Schmidt states. Or, we calculate the expectation 
value of an observable O-p, 

(6 V ) = TivhOv = Y, c p(p\Op\p) , (16) 



which is the expectation value of Op in a Schmidt state times the probability to find 
that state in the partonic subsystem summed over its whole basis. 

Thus, we conclude that via the Schmidt decomposition the subsystem and its 
environment can be described by simultaneously diagonal mixed state density ma- 
trices (/S-p g 7^ Pv,q)- In particular, note the absence of quantum interference terms 
oc (p\p)(p'\p) in eq. ([15]), for example, or oc (p'\Op\p) in eq. (|16|), which repre- 



sent incoherent sums. Hence, we achieved what is called quantum decoherence in the 
subsystem [0, |J. Note that in distinction to eq. (|T6| ) we obtain 

(6 g ) = Trgpgdg = £ c p (g\p n g ) (p n g \6 g \g) , (17) 

P,9 

where interference terms are present. In any case, the fact remains that the complex 
system is in the pure state |^): decoherence is achieved in the subsystem at the 
expense of establishing unique quantum correlations between subsystem and environ- 
ment which are explicitly constructed in eq. ( |Tl~D above. 

Next, focussing our interest on the partonic subsystem, we define its von Neumann 
or statistical entropy as usual . 



S-p = — Tr-p p-p In p-p = — c p In c p > , 



recalling that pp is diagonal with < c p < 1, see eqs. ( p4|14D , and assuming always 
at least two non-vanishing terms in the sum over states. Thus, a non-zero entropy 
emerges in a complex pure-state quantum system quite naturally, if only a subsystem 
is considered (or accessible to experiment) and the trace over the remaining environ- 
ment degrees of freedom is calculated to obtain the relevant density submatrix. Note 
that the total entropy vanishes, S = — Tr p In p = — 1 • In 1 = , as it should 
be. The same conclusions can be reached for the linear entropy employed in Refs. 



il 



s l 4 n = Tip {pp-pv) = i - tw># = i - Es 2 > ( 19 ) 
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which more directly measures the "(im)purity" of a density matrix (0 < Sip 1 < 1) and 



is often easier to handle than the standard statistical entropy. In Sec. 4.2, eq. (114), 
we show that the linear entropy provides a lower bound on the statistical entropy. 

Finally, the evolution of the complex system has to be taken into account. Up to 
now our considerations were limited to a fixed pure state and its Schmidt decompo- 
sition which, therefore, has to be recalculated from one instant to the next following 
the evolving system. Of course, the unitary evolution according to the Schrodinger 
equation, 

idt |#> = H |*) , |tf> = e"^ |#o> , (20) 

leaves the closed quantum system in a pure state, i.e. eq. (JIJ) remains valid with 
denoting the time-dependent state vector. We may split the hamiltonian, 

H = Hp + Hp g + Hg , (21) 

where Hp tg\ acts only on the partonic subsystem (gluonic environment) and Hpg 
describes the interaction connecting the two factorized subspaces introduced above, 
see eq. (Q). 

Until the end of this section we neglect the interaction term, Hpg « 0, and consider 
only an exactly separable system for simplicity. Furthermore, let us assume that a 
Schmidt basis of V simultaneously diagonalizes Hp. This implies 

[ p v , Hp ] = , and d t p v = . (22) 

Thus, we may write Hp = J2 p E p \p)(p\ , where E p are the associated energy 
eigenvalues. Then, the most general p-p allowed by eqs. ( p2|) and orthonormality of 
the Schmidt states is of the form given by eq. (]T2|). It includes the interesting case 
of a "pseudo-thermal ensemble" , 

Trp exv{-J2 P PpEp \P)\P\} p Z 

where Z(/3 P ) = J2 P ex V>{~PpE P }', for (3 P — > [3 the usual density operator of a canonical 
ensemble results [p3[| . Of course, we find a non-zero entropy here by eqs. ([T8| ) or 
(P^|). We remark that eqs. (P2"D by themselves imply that the expansion coefficients 
(probabilities) in eqs. (|l|,|13|) are constant, d t c p = 0. However, from our assumption 
that {\p)} diagonalizes Hp, we cannot further conclude [ pg , Hg ] = 0. In general, 
depending on the initial state |^ ) m e Q- (10)' the Schmidt basis states \pg) G Q will 
show a complicated time-dependence. 

Exact separability, naturally, cannot be expected to correspond to physically in- 
teresting situations in any realistic way: (A) Consider two decoupled field theories, 
one of which, the partonic subsystem, is additionally assumed to be essentially a free 
theory, cf. eqs. above; properly (anti) symmetrized partonic multiparticle states 



built from single parton momentum eigenstates constitute a suitable Schmidt basis 
here; coherent states do as well, if one neglects complications due to their overcom- 
pleteness. - Note that a Schmidt basis is, of course, not restricted to single-particle 
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states. - Or, (B) consider a metal as a free electron gas in a background lattice, which 



in turn can be described by a free phonon gas in the harmonic limit . 

At this point we anticipate that a non- vanishing interaction tl-pg will significantly 
influence the evolution of the subsystem and its environment. In Sec. 3 we present 
analytical results for a rather simple interacting complex system. It turns out that 
under favourable conditions, loosely speaking in terms of the previous example (A), 
the partonic Schmidt basis remains rather stable with respect to the evolution of 
the complex system and, especially, allows a parton to be represented by a slowly 
evolving gaussian wavepacket. Thus, it behaves essentially as a classical particle 
under the influence of the interaction with the gluonic environment. 

This result provides part of the motivation for our attempt to understand entropy 
production in strong interactions at high energy as being related to the quantum 
decoherence of the relevant Schmidt states, which would be stable in an idealized 
situation. Such states are commonly referred to as "pointer states" |7|, ^TJ . Typically 
a pointer in a classical measuring apparatus corresponds to a preferred quantum 
pointer basis being selected dynamically, such that the interactions in the complex 
system (observed subsystem + environment + measuring apparatus) preserve the 
stable classical character of the pointer. Its positions on a scale can be described by 
states, the quantum superpositions of which are dynamically suppressed (meaning 
decoherence). Thus, in distinction to the well-defined Schmidt states, pointer states 
refer to a dynamical situation. Up to now no general criteria exist which guarantee 
their existence. Sometimes they may be realized only in an approximate way in a 
given complex system [^TJ, which not necessarily has to be separable into a classically 
behaving subsystem plus its (quantum correlated) environment. The purpose of Sec. 
4 is to begin to study these questions for partons in a gluonic environment. 



3. Entropy Production via Pointer States in a Non-relativistic Model 



3.1 Time-evolution of the Density Submatrix 

In the following we study a complex system of a particle linearly coupled to 
oscillators as a simple example to illustrate the concepts of the previous Section 2. 
For a most relevant review of previous work on the Caldeira-Leggett model in the 
context of quantum Brownian motion we refer to Ref. PB| , see in particular Part III. 



Variants of this model have been applied before to the quantum decoherence problem 
("reduction of the wave packet") |24| and to the problem of radiation damping for 
a non-relativistic electron coupled to the quantized electromagnetic field [p7 |, where 



special attention has been paid to the initial condition of the complex system p8 |. 

For our purposes the Caldeira-Leggett hamiltonian can be written in analogy to 
eq. (|21|) with 

2 N ( 2 ~\ 

^-h' jr « s S{^: + * m ^} • (24) 

describing the partonic subsystem ( "parton" ) by a single free non-relativistic particle 
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(moving in one dimension for simplicity) and the gluonic environment by a set of N 
harmonic oscillators, and with the linear coupling given by 



H V g = -xJ2 C n\x n - J™ To. X ) ' ( 25 ) 
n=l I 



2m n ul 



Thus 



F = £ + £{£ + ^(*-s&*)'} ■ (26) 

Choosing c n = rrinU) 2 , we see that the model becomes explicitly translationally in- 
variant |26], |27], |28|], which presently also motivates the introduction of the parton 
"potential renormalization" term oc x 2 in eq. (|25[) ; it can be visualized as a "parton 
moving along with its gluonic springs attached 1 . In this form the model still allows to 
describe a large variety of physically interesting systems by choosing an appropriate 
spectral density for the environment []56 . 



N 2 i N 

As is well known, the environment induces noise and dissipation in the subsystem, 
which are essentially determined by I(uj) (and initial conditions). Assuming a quasi- 
continuous distribution of environmental oscillators, spectral densities I(oo) oc uj k (for 
sufficiently small u) with k > have been widely studied before, e.g. [|4], |26|, p8| . 
They correspond to so-called sub-Ohmic (k < 1), Ohmic (k = 1), and supra-Ohmic 
(k > 1) environments, respectively. - The electron-radiation field system involves a 
supra-Ohmic environment with k = 3 p7| . 



Here we consider a general spectral density in terms of an arbitrary dimensionless 
function F, 

I(lu) = g^F^u) e(fi - w) , (28) 

where g is a dimensionless coupling constant and Q denotes a high-frequency cutoff. 
For later purposes we introduce the (zero temperature) "noise" and "dissipation" 
kernels [24, 26], v and rj, 

roc /*oo 

u(s) = duj IUo) cos(o;s) , r/(s) = — / <iu; /(a;) sin (cjs) ; (29) 

JO JO 

see, for example, eqs. (BDlffTl) below for their dynamical effects. In the short-time 



limit, Qs -C 1, we obtain: 

= g u ^[l + O(fiV)] , V {s) = -g v Q*[Qs + 0{Q 3 s 3 )} , (30) 

with g u = g Jq 1 dx F(x) and g v = g Jq 1 dx xF(x). This limit may be completely irrele- 
vant for macroscopic quantum systems. However, assuming Q to be vaguely related 
to a separation of non-perturbative strong interactions with the gluonic environment 
from a perturbative regime at a scale on the order of Aqcd ~ 200 MeV and relevant 
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time scales considerably less than 1 fm, we expect that something analogous to the 
short-time limit will be of interest there. 

Having specified our toy model as above, our aim now is to calculate the density 
matrix p v of the partonic subsystem, see eq.(|), and to study the consequences of its 
time evolution. To achieve this, the dependence on the gluonic environment (harmonic 
oscillator degrees of freedom) has to be integrated out. In the Caldeira-Leggett model 
this can be done exactly using the Feynman- Vernon influence functional technique 



23 , |26| , since the coordinates (and momenta) appear at most quadratically in the 



environmental hamiltonian and in the interaction, cf. eqs. 

Assuming the initial condition that the total density matrix factorizes, 

p(0) = pp(0).p o (0) , (31) 

where p-p(O) describes the initial state of the parton, which will be defined below, 
and pg(0) denotes the density matrix of the gluonic environment which is assumed 
to be in the ground state at to — 0; i.e., the environmental oscillators only perform 
zero-point motion ( "vacuum fluctuations" ) initially. Then, the parton density matrix 
at a later time t is calculated with the appropriate propagator, 

p-p(x f ,x' f ,t) = J dxi dx[ J(xf,x' f ,t;xi,Xi,0) p-p(xi, x[,0) , (32) 



which has the path integral representation 26 



7(^,4,^,^,0) = LjvqVt/j&M-SrWe-^'fl , (33) 

with boundary conditions q(0) = x iy q{t) = Xf, and for q' analogously. Here Z is a 
normalization factor, which will conveniently be fixed at the end of the calculation, in 
order to preserve the normalization of p-p or, equivalently, of the parton wave function; 
S-p denotes the free parton action corresponding to Hp, eq. (p4|), 

S v [q] = fds\Mq 2 , (34) 



and $ is the Feynman- Vernon influence functional describing the influence of the 
environment on the subsystem. In the present case it is completely determined by 
the noise and dissipation kernels |2"6| , eqs. 

^[Q, q'] = [ ds [ du [q(s) - q'(s)] {u(s - u)[q{u) - q'{u)} + i-qis - u)[q{u) + q'(u)]} 
Jo Jo 

+ ^ioc[q,q'} , (35) 
where the "localizing" part of the influence functional, 

$ ioc [g, q'] EE i ^ \ ds [q 2 (s) - q' 2 (s)} , (36) 
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comes from the term oc x 2 in eq. ([25|). The constant appearing in eq. (|36|) can be 
calculated, 



- e 4j = I E ^ - f £ m - *f> • (37) 



where we used eqs. (^7|,^|) at intermediate steps and g = g Jq dx x~ 1 F(x). Equation 
([37] ) shows that I(ui) has to vanish sufficiently fast in the infrared, in order to avoid 
a divergence here. 

In any case, we observe that the resulting path integrals are gaussian and, hence, 
can be evaluated exactly. Introducing new coordinates (with a unit Jacobian), y = 
q — q' and z = |(g + q r ), one obtains from eq. (|33|): 

J(y f ,z f ,t;y h z h 0) = ^JVyVze 1 ^^^ , (38) 

with the boundary conditions y(0) = yi = Xi — x\, z(0) = Zi — |(xj + a^), and 
y(t) = yj — Xf — x'f, z{t) = Zf — \{xf + x'j), and where the relevant total action now 
is given by 

S[y, z\ = J dslftdyz — Muj^yz + J du y(s)[iv(s — u)y(u) — 2r](s — u)z(u)] 

(39) 

Since v is an even function, the second to last term here may be rewritten such that 
both integrations range from to t. 

The calculation of the path integrals proceeds in the standard way [^, |26[ . Since 
they are gaussian, the essential coordinate dependence of the propagator can be ob- 
tained from the total action, eq. (0), evaluated with the extremal ("classical") 
paths.0 They are described by the usual Euler-Lagrange equations. Separate varia- 
tions w.r.t. y and z, respectively, yield the equations of motion, 

M'z(s) + MunZ + 2 / du ij(s - u)z(u) = i I du u(s - u)y(u) , (40) 
Jo Jo 

My(s) + Mu%y - 2 J* du r](s - u)y(u) = . (41) 

Using the transformation y(t — s) = y{s), one finds that eq. (O) written in terms 
of y(s = t — s) formally coincides with the homogeneous part of eq. ([40]). Thus, 
starting at s — t, it describes a corresponding motion backwards in time, however, 
with different boundary conditions. Since the equations are linear, real and imaginary 
parts of the solution for z can be calculated separately (y is real by eq. flUf ) and the 
boundary conditions). Fortunately, it turns out that the imaginary part of z does not 
contribute to the minimal action [^, ^] and, therefore, is irrelevant at present. 

Then, we only have to solve one equation with the structure of the homogeneous 
real part of eq. (|40|) with appropriate boundary conditions. Its general solution can 



5 An overall time-dependent factor stemming from the remaining integrations over gaussian fluc- 
tuations around the classical paths is absorbed into the normalization factor Z, see eq. (p3). 
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be easily found in the short-time limit, cf. eqs. (30), with the Laplace transform 
technique, 

*00 = 1^-12 { z (°)lf- cosh(/_ S ) + f + oos(/ + «)] + z(0)if- srah(/_a) + /+ sm(f + s)]} 
J — + J + 

(42) 

where / ± = Q{±g Q/M + [(g n/M) 2 + 2g r) Q/M} 1 / 2 }^ 2 and g as introduced in eq. 
(|37|) . Employing the transformation mentioned after eq. fl4l|) and inserting the correct 
boundary conditions from eq. (05), we obtain: 

z(s) = Zi[g_ cosh(s^) + g + cos(s + )] 

, , s , M-, f-sinh(s_)+ f + sin(s + ) , . 

+ {„ - cosh( ( _) + 9+ co s(i+ )]} ^-A_L_^_^ , (43) 

= JZ/tfl'- cosh(i_ - s_) + cos(t + - s+)] 

/_ sinh(t_ — s_) + / + sin(t + — s + ) 



+ {!/* - cosh(t_) + cos(t+)]} 



/_ sinh(t_) + /+sin(t+) 

(44) 



with s± = /±s, t± = /±t and #± = f±{ft + Note that in the short-time 

weak-coupling limit < s± < t± € 1. In this limit the results could be further 
simplified. 

However, in the short-time limit generally we have to use the full expressions 
for the classical trajectories from eqs. ( f4~3| , f44"|) , when calculating the minimal action. 
As mentioned above, only the real part of z contributes to it. Thus, after a partial 
integration of the first term on the r.h.s. of eq. fl39|) and making use of the real part 
of eq. (flUl), the minimal action is given by 



Hmin(y/, Zf,t;yi,Zi,0) = M (yfZf - yik) + -i ds / du v{s - u)y(s)y(u) 



M (y f Zf - yik) + ^ig u Q 4 



ds y(s) 



(45) 



using (p0|) in the last step and where Zij = z(s = 0,t), as determined by eq. (§^). 
Inserting eqs. (H1,H3), the final result is: 



T, min (y f ,z f ,t;y u z h 0) = E R + iEj , (46) 
in terms of the real and imaginary parts of the minimal action, 

/_smh(t_H/ + sin(t + ) { lVfZf + VA] lf - COSHt - ] + fl C ° S(t+)] 

Vf^ld-f- + 9+fl + cosh(t_) cos(t + ) + (g+ - g~)f+f- sinh(t_) sin(t + )]} 

ay f Zf + by^i + cy f Zi + dy^j , (47) 
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— 2 {j/*[cosh(t_) - cos(*+)] 



[/_ sinh(t_) + /+ sin(t 

+yA^T7 ± sinh (*-) sin (*+) + (9+ ~ 9-)(l - cosh(t_) cos(t + ))]} 2 
J-J+ 

= Ay 2 + B yi y f + Cyj , (48) 

and where we introduced the abbreviations for the time-dependent coefficients for 
later convenience. We remark that for the non-interacting case, i.e. go = g v = g v = 0, 
the imaginary part vanishes, whereas the real part reproduces the well-known 1/t- 
term, 

s /? = —{yf z j + Vi z i-yf z i-yi z j} , ( 49 ) 

which characterizes the free particle density matrix propagator. Collecting the above 
results, eqs. (BS) and (HBT) - (ffH), the propagator J immediately becomes 



J{yf,Zf,t;yi,Zi,Q) = J7X^ R ~^ T • ( 50 ) 

The normalization factor Z (cf. the last footnote) will be calculated shortly. 

We consider, in particular, a normalized gaussian wave packet as the initial state 
of a parton with momentum p, 

ifj(xi,0) = (Xi,0\p) = ae~2 x i/^ 2 e ipXi , (51) 
with a corresponding density matrix, 

p P (xi,x-,0) = (Xi, 0^)^14,0) = a 2 e -\^ 2 i + x ?) 1 1 32 ~ x 'i) . (52) 

Then, normalization initially requires 

/CO 
dx p-p(x, x,0) = ir 1 ' 2 a 2 [3 , (53) 
-00 

which fixes a = 7r~ 1//4 /5 _1//2 . The normalization, however, has to be preserved through 
the time evolution of the system. Therefore, we proceed to calculate the time- 
dependent parton density matrix. It follows from eq. (^), using eqs. (|4T|,4*8,50,52|) 
and the coordinate transformation from eq. (^) above, after doing the two gaussian 
integrals, 

PT(yhZf , t) = *L .(^(W.'/.O/t)' + " (C+\PV) y} (M) 

with 

D(y f ,z f ,t) = idz f + ip-(B + l(3 2 bc)y f , £ = (A + \(3~ 2 + \(3 2 b 2 ) 112 . (55) 
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We remark that 

\i/;{x,t)\ 2 = f>p(x,x,t) = p v (y f = 0,z f = x,t) = t^-z-~M x + vfie ) (56) 
which, when integrated over x, fixes the normalization of the propagator. The result 



is: 



Z{t) = 2vr|d|- 1 , (57) 



which, of course, has to come out independently of the chosen initial state. 

Several interesting limits of the above results, which are generally valid in the 
short-time approximation, might be studied. It is instructive to evaluate the non- 



interacting case first of all. Inserting the time-dependent coefficients from eqs. ( f47| , f48 
we obtain from eqs. (]5"6[,57) in the limit g — > 0: 



£( yf = Q)Zf = X)t ) = n -l/2 w -l( t)e -(x-V t) 2 /w 2 Q (t) ; (5g) 

w (t) = (P 2 + P- 2 M-H 2 y/ 2 , * = £ , (59) 

which shows the well-known spreading of the wave packet due to the time-dependent 
width Wo and its motion with velocity Vq according to the classical law. By comparison 
with eqs. (^,^) we find the corresponding width and velocity for the interacting 

system, 

w(t) = 2£\d\- 1 , v(t) = -pd-H- 1 . (60) 
Thus, using eqs. ([54]) - (|57|) and (0), the full density matrix becomes 



Pv(y f ,z f ,t) = 7r-V2 w -i(t)e-(*/- v (*)*)V™ 2 (*) 

_ -y 2 (C + \{3 2 c 2 - d~ 2 [B + \(3 2 bcf/w 2 {t) 

e iy f (az f - 2d-\B + \(5 2 bc]{z f - v(t)t)/w 2 (t)) ^ (gl) 

which for g — > reduces to 

. jyfMt- 1 ( Zf [l - (3 2 /w 2 {t)] + v t(3 2 /w 2 (t)) ? (62) 



cf. eqs. (H| 

3.2 Quantum Decoherence in the Short-time Strong-coupling Limit 

Next, we turn to a discussion of the short-time strong- coupling limit, which we 
define through the following conditions: 

1. fit < 1 ("short time") ; 

2. g > 1 , g /g V) n > 1 ( "strong coupling" ) ; 

3. tt/M > 1 . (63) 
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Note that the second of conditions 2. implies that the function F specifying the 
spectral density J, eq. (^8|), has to be strongly infrared dominated by the definitions 
of g VtV and g following eqs. (|30| ) and (0), respectively. Condition 3. is added mainly 
for convenience. From conditions 1. - 3. we obtain 

fl « 2g (Q/M)Q 2 , /+»/_, « 1 , g- « 1 , (64) 



where we used the definitions of /± and following eqs. fl42| ) and (f44|), respectively. 
To be definite, we will be interested particularly in times t such that / + i < 0(1). 
Some approximate relations follow here for the time-dependent coefficients defined in 



eqs. A47L48 



= b « Mf + cot(f + t) , c « d w -Mf + / sin(f + t) , 
An C « iS « i^ 4 [l-cos(/ + t)] 2 /[/2 sin 2 (/ + t)] , (65) 



which will be useful henceforth. Note that all of eqs. ( pop approach the correct limit 
for g — > or t — > 0. Evaluating eqs. fl60D in the short-time strong-coupling limit, 
using (|65|) , one finds 

p sin(/+t) sin(/ + t) 

w m = Vo -j-r ' (66) 



w\t) « /3 2 [cos 2 (/ + t) + 2 2 ^ n)2 (l - cos(/ + t)) 2 ] + /3- 2 M- 2 /- 2 sin 2 (/ + t) . (67) 



Generally, the results show several remarkable features, if one compares the expression 
for the density matrix obtained from eq. ( |6~T1) in the present limit with the non- 
interacting case given in eq. (|62|). 

To begin with, all time-dependent functions entering eqs. (|60| , |6~TD through the 
coefficients evaluated in eqs. ( 165] ) are governed by a single dynamical time scale f^ 1 , 
cf. (|64]). From previous experience we expect that their precise form is sensitive to 
the initial condition [21, [2"8[| , particularly the simple factorization in eq. (|3"T|). The 



fact that in the present limit trigonometric functions dominate over the hyperbolic 
ones in eqs. (|47l ,|48|) , which usually would rather determine the long-time behaviour 
of the propagator etc., also points towards this conclusion. We remark that for t — » 
eq. (|6l|) yields p-p — > pp, eq. (|62"D, in agreement with the factorized initial condition. 
Furthermore, there is obviously and quite naturally a dependence on the initial state 
of the parton parametrized by the width j3 and momentum p. Note, however, had we 
chosen, for example, a gaussian superposition in momenta instead of eq. (|5"T|), 

0) = a e-^/P 2 [°° dp' jP'^ e~(P' ~ ^ I p2 , (68) 



this would simply replace £ 2 — > ^ 2 + |P 2 , cf. (|55|), and leave the density matrix 
the same otherwise, eq. (p4|) , apart from a constant overall factor. We will consider 
another, more interesting parton initial state below. 
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Turning to the diagonal elements of the density matrix, p-p(0, Zf, t) from eq. ( |56"D 
and (^l|) in particular, we observe that they are completely determined by the velocity 
and width, v and w as obtained in eqs. ( |66"D and (|6T|), respectively: 

(i) The free velocity t>o is modified by a slowly decreasing function of time (scaled by 
/+), which shows an influence of the gluonic environment oscillators on the parton as 
an effective friction force. 

(ii) Compared to the spreading of the wave packet of a free particle, see wq in (|59|) , the 
width has become a very slowly evolving function of time in the region of interest to 
us, f + t < 0(1), and can actually be smaller than the initial width w(0) = wo(0) = f3 
for a suitable choice of the parameters {localization}] note that the term oc g u , which 



is due to the noise kernel, cf. (|29"D, presents an increasing contribution to w which 
can be neglected, if the initial width is not excessively too small. 

In fig. 1 we illustrate the above results by showing the unique deceleration func- 
tion v/v , see eq. (|66"D, as a function of the dimensionless variable t + = f + t = 
ilt(2g Q Q,/M) 1 ^ 2 . Furthermore, we show the relative change of the width of the parton 
wave packet, cf. eq. (|67|) , as a function of the same variable: 

- 13 = {cos 2 (t + ) + a<[l-cos(t + )] 2 + asin 2 (t + )} 1/2 -l , (69) 

with the parameters given by 

a = ^T,W > «< = W • (70) 

2fi- M g il 



Note that the expression corresponding to eq. (|69D for a free particle (g = 0) can be 
written similarly, 

^iMzJL = |i + at+ 2} 1/2 _ 1 , (71) 

which is also shown in one case in fig. 1 for comparison. 

By modifying the spectral function of our toy model, eq. (|28D , one could opti- 
mize the above results in the sense that (i) the environmental friction is minimized 
and simultaneously (ii) the stability of the parton wave packet is maximized. It is, 
however, a general property of linear coupling models a la Caldeira-Leggett that fric- 



tion and localization go together to some extent ||26|| . We remind ourselves that for 
supra-Ohmic environments with spectral exponent k > 2, as discussed after eq. (fZ7j), 
one even finds an effectively free quantum particle with renormalized properties in 
the long-time limit. Interestingly, this applies for a non-relativistic electron coupled 
to the quantized electromagnetic field, where k = 3 |2jJ . Thus, as could be expected, 
our parton toy model specified by eqs. (p8|) - ( |30|) and (0), in particular, behaves 
very different from the electron-radiation field system. 

At this point we want to pause for a moment and turn our attention back to the 
discussion of partonic Schmidt and pointer states in Sec. 2 before. After all, the toy 
model in this section is meant to illustrate those concepts and to provide an explicit 
example, how they may be realized in a physical system, which is vaguely oriented to- 
wards QCD partons. It should have become clear that the single-parton wave packet, 
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eq. (|5lD, trivially constitutes a Schmidt basis for the single-parton (plus gluonic en- 
vironment) system under consideration, cf. eqs. ( [TT]) - fll3]). However, far less trivial 
seems the fact that this Schmidt basis also is dynamically rather stable in the sense 
of points (i) and (ii) above. Thus, a parton here behaves in a good approximation 
as a classical particle moving in a dissipative environment. Hence, we have shown 
that the Schmidt basis defined by the chosen initial wave packet state also presents 
an approximate realization of a pointer state basis. To complete this demonstration, 
we have to show that these pointer states have appropriate decoherence properties as 
discussed after eq. (|16|). 

Thus, instead of considering the off-diagonal density matrix elements from eq. 
(|6l|) (for Uf ^ 0) in detail,^ we now consider a different parton initial state. Our aim 
is to show that 

(iii) the gluonic environment enforces the evolution of a pure state consisting of a 
coherent superposition of pointer states into a mixed state consisting of a decoherent 
sum related to those states; i.e., the interference terms of the initial density matrix 
become dynamically suppressed as time goes on. 

This is the origin of entropy production in the present model, as we shall see. 
Instead of the parton initial state defined in eq. (|5"ID, we want to study the 



coherent superposition of two such wave packets with momenta +p and — p: 

fe.O) = iV- 1 / 2 [^(x l ,0) + V-p(x 4 ,0)] , N = 2[l + e-(P0) 2 ] , (72) 

where the index ±p is introduced to indicate the respective momenta and N is the 
necessary additional normalization constant, see eqs. ([51]) - ( [5"3|) for comparison. 
Equation (|72"|) yields the corresponding initial density matrix, 

= N~ 1 [p v (x i ,x' i ,0)+p + p v (x i ,x' i ,0)_ i) + pi„i(xi,x-,0)] , (73) 

with the pure state density matrix elements p-p(xi, x[, Q)±p as in eq. (0) and the 
interference term 

Pnizi^O) = Q J e-K a ?+ a ?)//' a [e* + !C <) + 6* + ^] . (74) 



Then, the time evolution of p-p, which follows the linear law of eq. (p2|) , can be 
evaluated term by term. The result for p v ±p is given by eq. ( pTI ) with the appropriate 
momentum inserted, whereas the interference term at a later time t is calculated to 
be 

p int (y f ,z f ,t) = Y^Pv{y f ,z f ,t)\ v ^^ cyf ^ c (y f±2 p/c) ■ ( 75 ) 

Here p-p is given by eq. (EHl) as well, however, with the indicated substitutions, which 
can easily be understood. Since presently v — *■ 0, i.e. p — 0, there is no motion 

6 Generally, one finds with the help of eqs. (^) an increased gaussian (in yj) fall-off as compared 
to the non-interacting case, eq. (|62|), with time-dependences governed by the scale /+ 1 discussed 
above, and corresponding modifications of the phase factor. 
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of the gaussian in the center coordinate Zf, compare eq. (|61|). On the other hand, 
the substitutions q// — > c(yf ± 2p/c) corresponding to the respective parts of the 
interference term, eq. (|74[), imply that parts of the resulting gaussian in the relative 
coordinate yj, see eq. flSTD , "move away" with a velocity proportional to the relative 
momentum 2p of the superposed initial wave packets. There is also a similar effect 



on the phase factor. Because of the complexity of the r.h.s. of eq. (|75), it is not 
illuminating to examine the behaviour of individual density matrix elements obtained 
from eq. (|73[) at a later time t with the help of eqs. (|6TD and (fT5|). 

Therefore, we proceed to calculate directly the entropy for the various interesting 
states. In particular, we employ the linear entropy defined in eq. for two reasons: 
a) S!p n has a simple bound from above, which eventually gives a good qualitative im- 
pression of the saturation of entropy production; b) practically, for any somewhat 
complicated non-diagonal density matrix, it seems impossible to calculate the stan- 
dard entropy S-p according to eq. (|IS|) without performing a more or less difficult 
diagonalization first of all. 

Let us begin by calculating the linear entropy for the simple pointer state 



(pl|) with its density matrix ( |61~D using eq. (|T9|): 

/oo /*oo 
dy f / dz f p v {y f) z f , t) Pr(-y f , z f , t) 
-oo J —oo 

= l-^c^w- 1 , (76) 

Cl = C + \(3 2 c 2 - d- 2 w- 2 [B + \f3 2 bcf , (77) 

independently of the parton momentum p. In the non-interacting limit (g — > 0), 
where the previously defined time-dependent functions reduce to A = B = C = 0, 
a = b = — c = —d = M/t, and w wo, v — > Vo, eqs. (|59|), the r.h.s. of eq. ( [76]) 
correctly goes to zero for all times. This is in agreement with eq. Ql9D , since here the 
initial pure state remains a pure state, which has zero entropy. However, our proposed 
pointer state (|5l|) , is not perfectly stable. Consequently, we expect some amount of 
entropy to be generated due to the interaction with the gluonic environment. Indeed, 
in the short-time strong-coupling limit (|63|) using eqs. (|65|) , one obtains from eq. ([76]) 



-> , < , i i 2 



ohn(j.\ ^ iy v 



2^ [l-cos(/ +t )P + ^si„=(/ +t) 



(78) 



which is very small for suitable f3Vt = 0(1) and according to the assumptions flB"5|). In 
this sense the pointer state (O) conforms with the so-called predictability sieve [B3 



which characterizes pointer states as producing the least possible entropy increase. 
Furthermore, note that a perfect pointer state would have to diagonalize the relevant 
part of the total hamiltonian, H v + H-pg, as defined in eqs. (|2^j25|) . We did not 



seriously try to improve the wave function ( |5T| ) to optimize its pointer state character. 
However, apart from the points (i) and (ii) discussed above, we will find additional 
support for having made a reasonable guess by its decoherence properties. 

For this purpose we calculate the linear entropy S ™ for the coherent superposition 
of two pointer states, eq. (]72|), with its density matrix p at a finite time t following 
from eqs. (|61]) and (p^) - (|75|) as described after eq. ([741). The result of a lengthy but 
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straightforward calculation, which involves a number of double gaussian integrations 
as in eq. ( [76|) , can be represented as follows: 

S l i n = 1-Tr v p 2 = l-N~ 2 Ti v {p v+ p + p P ^ + p mt } 2 = (79) 
1 _ iV-V 172 ^ 1 {[2-| + e - 2 ^AM 2 - \ci/ Cl] + [e 2c 3 + |c 2 2 / Cl + e -2(^) 2 ] 
+ [4 cos(|c 5 p/rf + ic 2 c 6 / Cl ) e -i^/^) 2 + & + 8c 3 c i " c e 2 " ^ 2 c 5 2 Ci)/8 Cl 



where we grouped terms into square brackets which come from (p-p + p + p-p -p) 2 , p 2 nV 



and the cross terms involving p int at last, respectively; N was defined in eq. (|T2|), C\ 
in eq. (|77|) , and the other abbreviations are 

c 2 = -cp(3 2 [l - (3 2 b 2 d~ 2 w- 2 -2Bbc- 1 d- 2 w- 2 } , c 3 = -p 2 /3 2 [l - /^Vd" 2 ^ 2 ] , 
c 5 = -2pf3 2 bd- l w- 2 , c 6 = -2pd~ 2 w~ 2 [B + \l3 2 bc] . (80) 

We should, of course, check again that in the non-interacting limit there is no entropy 
produced, since the pure state tp, eq. ([72]), which corresponds to p, remains a pure 
state. In fact, it is instructive to see, how this result comes about, 

S l 2 % = l-2iV- 2 |[2-i + e- 2 (^) 2 ] + [l + e- 2 (^) 2 ] + [4e-(^) 2 ]} = , 

(81) 

where the terms in square brackets are in one-to-one correspondence with those in eq. 
(|79|). We observe that all terms are exponentially small except for three contributions: 
two identical terms in the first bracket, which can be interpreted as a decoherent sum 
for the two superposed pointer states, cf. eqs. ( |76i[79l) ; however, there is also a 
large term originating from the square of the interference term, p 2 nt in eq. ([79|), 
which appears in the second bracket. Obviously, if there is any considerable entropy 
production to occur in the interacting case, then particularly this contribution of the 
interference term has to be significantly reduced, which essentially amounts to the 
decoherence effect that we are looking for. 

Using the full result, eqs. (^,77), evaluated in the short-time strong-coupling 
limit and expanding consistently to lowest non-trivial order in quadratically small 
quantities, i.e. terms which are proportional to gv/g 2 . or g u /(g Q/M), we find: 

S l 2 in (t) w S[ in (t) - 2N~ 2 ((Tx e ~ 2 ^) 2 - a 2 + 4a 3 e - ^)') , (82) 



where we also used eqs. ([76]) - (|78| ) and assumed f3Vt = 0(1) again; furthermore, the 
terms oc cri,2,3 correspond to the square brackets in eqs. (|79|) or (^l|) in an obvious 
way, with 

a, = (d-a^/A . (83) 
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Several remarks are in order here. Firstly, note that all cr's are oc p 2 . Thus, p — > 
correctly implies (t) — > S l i n (t), since in this limit we recover the simple (imper- 
fect) pointer state ([51]) with p = and its associated entropy, see eqs. ( [TED - ([78]). 
Secondly, for moderately large relative momentum between the two coherently super- 
posed pointer states in eq. ([72]) as compared to their spatial width, i.e. p > we 



obtain from eqs. (82,83) the final result 



S% n {t) w |a 2 . (84) 

Thus, we conclude that an appreciable amount of linear entropy can be produced here 
through the mechanism of environment-induced decoherence of coherently superposed 
( approximate ) pointer states. In particular, note that the dominant contribution oc o~2 
comes from the decay of the interference term, which we anticipated in the discussion 
following eq. (HI]). The rate is determined by the dynamical time scale cf. (|64T), 
which also governs the friction and localization effects mentioned in points (i) and 
(ii) above. 

It is worth noting that in the limit of very large relative momentum, p/3 —>■ oo, we 
find from eq. (|75|) in the short-time strong-coupling limit instead of fl8"l|): 

5f(t) « l-\[l-S^{t)][2.\] = i[l + • (85) 

This result corresponds simply to the expected decoherent superposition of two equally 
probable states, cf. eq. flT9|), which would yield S!p n = 1 — {(|) 2 + (|) 2 }- Again 
there is a small correction due to the imperfection of our pointer states as discussed 
above. 

Finally, we remark that for the demonstration of entropy production or environment- 
induced decohence we did not necessarily choose the optimal superposition of approx- 
imate pointer states in the sense of the largest or fastest effect. Our above result, 
eqs. flS2l) - fl84]) , shows that the entropy production is larger for superposed states 
which are further separated in momentum, i.e. states the wave packets of which have 
consequently moved further apart in space after a given amount of time has passed. 
From this observation and the fact that the interaction between partonic subsystem 
and gluonic environment in our toy model, see eqs. (^,^), acts on the respective 
coordinates rather than the momenta one may be tempted to speculate that the su- 
perposition of initially spatially (rather than in momentum) separated pointer states 
would decohere even faster and, thus, produce entropy at a higher rate. Indeed, in 



the different context of Ref. p4]j , such an effect of a preferred observable (coordinate 
vs. momentum) has been found. 

These remarks complete the demonstration of dynamically enforced decoherence 
of coherent superpositions of pointer states, point (iii) above, and its relation to 
entropy production for our toy model. 

More general superpositions of states can, of course, be studied similarly. In 
particular, a general parton initial state could be decomposed into coherent sums 
over pairs of localized gaussian wave packets (see also the discussion in Ref. for 
which the present considerations apply again. 
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The extension to a multi-parton system is straightforward only as long as we 
can neglect perturbative interactions among themselvesf] and interactions mediated 
through the self-coupling of the gluonic environment. Otherwise, provided that a 
Feynman- Vernon type influence functional can still be calculated (cf. Sec. 3.1), one 
has to face at least a generally quite difficult mechanical n-body problem, i.e. a large 
set of coupled equations replacing ( |40||41~1) . This seems to prohibit a stepwise devel- 
opment of our toy model into a more realistic phenomenological model. Therefore, 
we propose a fully relativistic quantum field theory model instead in the following 
section, which is motivated by the non-relativistic model studied here. It might mimic 
some features of QCD and the parton model. 

Generally, we expect to encounter the set of characteristic formal problems posed 
by the requirements of Lorentz and local gauge covariance as well as by the occurence 
of the usual divergences of quantum field theory. A more subtle difference, which is 
due to going over from the non-relativistic Schrodinger equation to proper wave equa- 
tions describing the basic fields, will mostly be hidden in our functional Schrodinger 
picture approach. 



4. Decoherence and Entropy Production in a Scalar "Parton" Field Theory 

4.1 Model of Momentum Space Mode Separation 

Before attempting to generalize our previous considerations to QCD gauge theory, 
we want to take an intermediate step here by studying a much simpler model of two 
coupled real scalar fields, 4>i an d 02, representing partons and "unobservable" gluonic 
environment modes, respectively. It is defined by the classical action 

S = J d 4 x {Li + L12 + L 2 } , (86) 
with (g^ = diagfl, -1, -1, -1]) 

Lj = |(%) 2 -^(0,) , J = 1,2 ; L 12 = -VfaM , (87) 

i.e., we consider an attractive interaction of the two scalar fields, in particular. Clearly, 
the meaning of the model depends on how we define the separation of partons and 
environment modes. 

To begin with, let us assume that there was only one scalar quantum field ip 
representing the total system. Then, similar to the background field method, we may 
split it, cp ee (f>i -\- (p2, according to the following considerations. 

The essential feature of QCD to be incorporated here phenomenologically is the 
running coupling constant. "Soft" interactions with small four-momentum transfer 
(squared, \Q\ 2 ) have a strong effective coupling and, generally, cannot be described 
in perturbation theory. Conversely, "hard" processes with large \Q\ 2 are reliably 

7 This is a question of the relative sizes of a non-perturbative time scale corresponding to 
here and a perturbative "rescattering" time, which has to be taken up in the QCD context 
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accounted for by perturbation theory and constitute the only basis for the comparison 
of QCD theory with experiment [ 29| . Field modes which are sufficiently close to 



constant field configurations and would interact only among themselves form a sector 
of the total Hilbert space which would be completely unobservable. They can neither 
hadronize nor can any large-|Q| 2 process occur within this sector due to the restriction 
of the four-momenta. Roughly, we assume fc 2 < m T 2 and k 2 < A.q CD . Of course, 
there is no sharp boundary of this sort in QCD. In any case, this defines the gluonic 
environment described by the field 02 and L2 in eqs. (|36|,f37]) and the region II in fig. 
2. (We don't distinguish soft quark modes presently.) The complement of this part of 
the total Hilbert space, region I in fig. 2, is described by 0i and L x and constitutes the 
parton subsystem, which is almost directly observable by hadronization or via large- 
\Q\ 2 processes. Finally, we describe the interactions between partons and environment 
by L 12 in eqs. (gEp). 

Several qualifying remarks are in order here: 

1. The Momentum Space Mode Separation ("MSMS") as introduced above is not 
Lorentz invariant. Therefore, it only makes sense in a certain frame. It should be 
the one in which we formulate the initial conditions for the time-evolution problem 
of our field theory (cf. Sec. 4.4). 

2. In a non-abelian gauge theory MSMS is not gauge invariant. Therefore, it only 
makes sense, if the gauge is completely fixed; particularly, if the theory is described 
in terms of physical degrees of freedom. 

3. At present 0i and 2 are quantum fields, whereas usually in the background field 
method the background consists of a classical field, which does not propagate in 
quantum loops, whereas the quantum field only lives on internal lines in Feynman 
diagrams by construction. Obviously, a separation of field modes according to some 
dynamical criterion (e.g. MSMS) requires a non-trivial modification of the standard 
background (gauge) field method [^(J . 

4. There is an analogy to the non-relativistic model of Sec. 3; the interaction L12 
may be depicted as the "partons dragging along the unobservable background field 
modes" or vice versa. A close analogy would be a non-local interaction involving a 
form factor, 

d A x L12 = J d A x d 4 x' [0i (x) — 02 (a;)] lu 2 (x — x') [0i(a/) — 02 , 

which is, however, local in momentum space and, thus, would rather directly gener- 
alize the interaction in eq. (p6j) . Unfortunately, such an effective interaction presents 
difficulties, when attempting the usual equal-time quantization of the model. There- 
fore, it is also not obvious how to write down a functional Schrodinger equation for 
it, cf. eq. (^) below, and we will not consider it any further. 

In order to specify our model, we consider the potentials in eqs. (^) in more 
detail according to MSMS. Renormalizability and stability require the potentials to 
be fourth order polynomials. Then, a minimal version of the model may be defined 
by: 

^l(0l) = ^| A 1^1 4 > ^(02) = -^M 2 02 2 + J| A 20 2 4 , 
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Vtfufa) = \u<f>?<f>i , (88) 

with the restriction < Ai < A 12 < A2, which is supposed to mimic in a very 
schematic way the | (^-dependence of the strong coupling. In t> 2 (0 2 ) we included a 
negative "mass" term to allow for a condensation of environment modes (generating 
different masses for <pi and 2 )- The 2 -potential is fully consistent with MSMS. 

However, in the mutual coupling we encounter the problem that it admits, for 
example, the annihilation of two space-like partons into two environment gluons with 
arbitrarily small \Q\ 2 - This contradicts a clear-cut MSMS, such that the respective 
size of couplings in ( [38l) corresponds to the |Q| 2 -dependent situation in QCD. The 
same effect would arise in other couplings, which we presently did not include for 
simplicity, and cannot be avoided in such a crude model. Similarly, the ^/-interaction 
admits some low-|Q| 2 processes, which it should not. Therefore, only on the average 
(in multiple scattering) may we hope to mimic a bit of the QCD case. This will be 
studied in the sequel to the present work. There we may want to include also 0i0 2 2 ~ 
and 0^02-couplings. They correspond to decay/fusion and soft scattering processes, 
which are important in parton cascades flTTl , pL2|, |18| . In the time-dependent Hartree- 
Fock approximation considered below the tree level Feynman diagrams for the two- 
point functions corresponding to the latter cubic couplings would simply carry one 
classical field (^1,2) on one leg, whereas l^(</>i, 2 ) as in (|88|) describes effective mass 
insertions (0j 2 )j=i,2 and something like a formfactor (</>i</>2), which we discussed above 
and which we consider most interesting to start with. 

Generally, the self-interactions Vj(<fij) and the mutual coupling V(4>i,(j)2) necessi- 
tate a renormalization of our model. It is well known, how to carry out the standard 
program for scalar fields. Furthermore, it has been extended to the study of initial- 
value problems recently |3tJ, which would be relevant in the following. However, 



since it does not illuminate our main topic here, we simply assume for the moment 
a regularization of the arising divergences (cf. below) by a high-momentum cut-off. 
In this way we also avoid a potential triviality problem (in four dimensions) of the 
quartic interactions. 

To proceed in analogy with Sec. 3, we consider the time- evolution problem for 
the quantum field theory defined by eqs. (|86|j87|) starting with given initial field con- 



figurations. In general, this can only be done in some approximation, e.g. ordinary 
perturbation theory for sufficiently small couplings in connection with any formal- 
ism able to handle time- dependent (non-equilibrium) systems such as (semi-classical) 
transport theory J3T[] or the Schwinger-Keldysh formalism (see e.g. Ref. |32| and 
references therein). Since the interesting results of Sec. 3 depend crucially on the 
short-time strong- coupling limit, cf. Sec. 3.2 and (^), in particular, and we an- 
ticipate the strong coupling at small four-momentum transfer in QCD to be most 
important for the environment-induced decoherence effects, we need a calculational 
scheme which allows us to consider the corresponding limit here. 

4.2 Variational Approach to Time-evolution in Quantum Field Theory 

It seems most efficient to employ the time-dependent variational principle of Dirac 



15, 16] which is equivalent to the functional Schrodinger equation describing the full 
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dynamics of a generic field (p in the Schrodinger picture, 



m\<P\t\ = H[%M*W\ = j d d ^{-\j^i + |(v^) 2 + v{<p)}v[r t t] , (89) 

where ^[^t] = ((p\^f(t)) denotes the wave functional in the (^-representation, which 
corresponds to a wave function ip(x,t) = (x\ip(t)) for a one-dimensional quantum- 
mechanical system, and tt = —i5/5(p represents the canonical momentum opera- 
tor conjugate to the field ("coordinate") (p. The dynamics is determined by the 
hamiltonian H. In this context the completeness and inner product relation, respec- 
tively, involve functional integrals instead of ordinary ones (orthogonality needs a 
^-functional), 

(*!(t)|* 2 (t)> = fV(p{^{t)\(p){^ 2 {t)) = JvtpVWwt^W] , (90) 
which induces 

(*(t)|0(7h jV »)|*(t)> = jvcp**[cp;t)0(-i—,cp)*[cp;t) , (91) 

i.e. the evaluation of expectation values of functions of operators corresponding to 
the usual quantum-mechanical formula 

Furthermore, note that eq. can be rewritten in terms of a density functional 

operator p, 



(¥(t)|O(7r,0)|¥(t)) = JV<p{<p\*(t))(*{t)\0(*,<p)\<d = Trp(t)0(7r,^) , (92) 
which is again analogous to the usual result, cf. eqs. (|T].[TB],[T7|) , for example. Finally, 



we may state the variational principle [30, 33 

^ffl = , for all with (tf(t)|tf(t)) = 1 , (93) 

and r[tf] ee [dt(V(t)\[idt-H]\ty(t)) , (94) 



i.e. requiring the effective action T defined in eq. Q9"l|) to be stationary against ar- 
bitrary variations of the normalized wave functional which vanish at t — > ±oo, is 
equivalent to the exact functional Schrodinger equation, eq. ( |59"D above. With the 
variational principle in hand, eqs. ( ]93|j9^ ), one can solve the time-evolution problem 
in quantum field theory approximately by restricting the variation of the wave func- 
tional \1/ to a subspace of the full Hilbert space, which means by choosing a suitably 
parametrized trial wave functional. 
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In the following we choose to work with most general gaussian trial wave junc- 
tionals. For the generic field if it is defined by 

Vafot] = N(t) expj-^-^^Gf-^-iSWjb-^)] + i*{t)[<p - <p(t)]} , 

(95) 

where here and henceforth we use a shorthand notation suppressing integrations over 
spatial variables, whenever they are obvious. For example, 

ypG-^t)^) = Jal d x d d y V {x)G-\x,y,t)${y,t) , 

7f(t) <p = J d d x 7t(x,t) tp(x) , and tr E(t) = J d d x E(x, x, t) .(96) 

The normalization factor iV can be easily calculated (for symmetric and positive- 
definite G) according to eq. (j9"0D, 



N(t) = (ATdet G(t))- 1/4 , 



(97) 



which is a useful result for further calculations employing (A/ - is an infinite constant 
which can be omitted in the following). The meaning of the variational parameter 
functions <f, ft, G, and £ follows from |3(], |33|j : 



(<p(x)) G = <p(x,t) 
(<p(x)(p{y)) G 



(-<■ 



5 



) G = ir{x,t) , 



5(f(x 

<p(x,t)<p(y,t) + G(x,y,t) 



(id t ) G = n(t)0(t)-ti[±(t)G(t)] 



(98) 

(99) 
(100) 



where the operator expectation values (■ ■ -)g are calculated according to eq. J9~lp with 



^g- (There is a trivial error common to Refs. |]30| , |33| in the sign of the term oc tt in 
the definition of the trial wave functional.) Thus, the trial wave functional *&g, eq. 
(p5|), is a gaussian centered at ft with a width G; ft and S are "conjugate momenta" 
for ip and G 8 , respectively. We anticipate that the application of a gaussian trial 
wave functional amounts to a description of the field theory under consideration in 
terms of coupled equations for one- and two-point Wightman functions ("TDHF", 
time-dependent Hartree-Fock approximation). The equivalence with the Cornwall- 
Jackiw-Tomboulis generating functional (effective action) for two-particle irreducible 
graphs was demonstrated by those authors in Ref. for energy eigenstates of the 
field. 

Next, we evaluate the effective action, eq. with the trial wave functional, eq. 
(1951), to obtain: 



r[* G ] = jdt{[v(p-\^-\(yvY-v{f)] 



8 This is obvious after a partial integration of the last term on the r.h.s. of eq. (100), which can 
be performed in the evaluation of the effective action, see eq. (101), for example. 
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+ tr 



EG - 2EGE - -G- 1 + -V 2 G 
8 2 



iy>;rV (2 ^)G(x.x) 
"l! V(4) ^ / ^^ 2 (x,x)} , (101) 



where we suppressed spatial integrations as in (|96|) where possible and the time- 
dependence everywhere, V 2 acts only on either one of the two arguments of G, 
V {n) (p) = d n V((p)/d<p n , and we assume the potential to be at most quartic. Note 
that the calculation of (V(<p))g, m particular, proceeds conveniently through a Tay- 
lor expansion at an intermediate step: 

<* O |V(¥0|¥ o ) = N 2 Jvp {v(^) + ^V( 2 )(^ 2 + iv( 4 )(^ 4 }e-^ G "V , 

(102) 

where we used eqs. (OT,pH) and observe that only even powers in (p contribute. 



Employing eq. (|97p , one obtains the simple gaussian integrals needed here by a 
suitable functional differentiation w.r.t. G _1 (cf. also eqs. ( 120 , 12l|) below), 



(p m (x)) G = N 2 jvipip 2 W{x)e-\vG V = i[3].GW{x,x) , (103) 



and, thus, produces the potential contributions to the above result, eq. ( |101|) . We 



remark that the terms in the first line on the r.h.s. of eq. (|101|) present the usual 



classical action, cf. the hamiltonian in eq. (jS9|) , whereas the terms in the second and 
last line are 0(h) and 0(h 2 ) quantum corrections, respectively. From the effective ac- 
tion the relevant coupled equations of motion are obtained by independent variations 
w.r.t. to the variational parameter functions <p, 7f, G, and E. This will be performed 
for our model in the following. 

4.3 Model-independent Entropy in TDHF Approximation 

After the above short review and formal preparations we are finally ready to study 
the scalar parton field theory defined in eqs. flSfiHTTj) with the variational technique. 
We assume a most general gaussian trial wave functional in product form, 

*i 2 [0i,<M = N 12 (t) ^ Gl [<f>i;t] VgAM 

■ exp {-§[0i- Mt)] [G 12 (t) - zE 12 (t)] [0 2 -fc(t)]} , (104) 



with the normalized gaussians on the r.h.s. as defined in eq. (p5|). Here, however, each 
of the variational parameter functions in the expression for *f>Q.[(f)j] t] carries an index 
j and A^ 12 denotes an additional normalization factor. The latter is necessary, since 
we included here an essential exponential describing possible two-point correlations 
between the parton field <pi and the gluonic environment field 02- It is obtained from 
the normalization condition, cf. eq. (j97[ ), 



1 = JvfaVfcVUfafat^ulfafat] 
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= {N x N 2 N X2 f J V^Vfo exp{-i[0 1 GTVi + 02{G' 2 - 1 -G' 12 G' 1 G' 12 }0 2 ]} 

= (N x N 2 Ni 2 ) 2 (det G^y 1 ' 2 (&et{G 2 l - G^dG^y 112 
= N,l (det{l - G 12 G 1 G 12 G2}y 1/2 

— > N 12 (t) = (det{l-G 1 (t)G 12 (OG' 2 (t)G 12 (t)}) 1/4 , (105) 

where we suitably shifted the fields to reach the second equation, assume G Xt2 to 
be symmetric and positive definite and similarly have to restrict the argument of 
the determinant in the final result to assure the existence of the functional integrals 
(dropping irrelevant constants as mentioned after eq. (|97|)). We remark that the 
"1" in eq. ( |105| ) represents a ^-function (of the space coordinates), which will occur 
frequently in the following. Furthermore, note that N X2 is symmetric under G± <-> G 2 , 
as it should be. 

To appreciate the importance of the non-diagonal factor (in <pi )2 ) in eq. ( |104| ), 
which involves parton-environment correlations, we calculate the partonic density 
functional according to the general definition discussed in Sec. 2, cf. eq. (^): 

pr(t) = Tr g p(t) = TV 2 |*i 2 (t))0Mt)| , (106) 

which yields the "matrix elements" 

(^IpviWi) = Prtthfat] = [vfcVl^fat^uttufat] (107) 



N&t) n^'iA^G^i-A 

Vfo V* G2 1<P2; t]^ G2 [0 2 ; t] exp 0i; iftfa - 2 )} 

= N^^M^gAM exp{|X 1 [0i,0 1 ;t]G 2 (t)X 1 [0i,0 1 ;i 
with 

X^^-t] = i[0i-0 1 (t)][G 12 (O+iS 1 2(t)] + |[^i-0i(*)][Gi 3 (O-iS 12 (t)] , 

(108) 

and where we used eq. (|104j) , see also eqs. ( p2| , |95|) . In addition to G X2 , cf. eq. 



( |105| ), we now also assume £i 2 to be symmetric, such that X x l = X x w.r.t. space 
coordinates. Obviously, if the correlation functions G X2 (t) and Si 2 (t) vanish, then 
the partonic density functional p-p is just the one of an independent scalar field in the 
pure state \^Gi(t))- According to our previous general considerations in Sec. 2, cf. 
also the application in Sec. 3.2, the partonic entropy has to vanish in this case. 

We can verify this conclusion by a direct calculation (similar to the one above) in 
the case of the linear entropy, which was defined in eq. (fL9|), using the parton density 
functional from eqs. ( |107| , |108] ): 

S l p = 1 - Tr r p£ = 1 - N& det(A)- 1 / 2 det (A - BA^B)- 1 ' 2 , (109) 
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where again, as in N% 2 from eq. (|105|) , we have to assume that the arguments of the 



determinants are symmetric and positive definite. They are defined in terms of 

A{t) = l-|G 1 (t)G 12 (t)G 2 (t)G 12 (t) + |G 1 (t)S 12 (t)G 2 (t)S 12 (t) , 

B(t) ee \G x {t)G 12 {t)G 2 {t)G 12 {t) + piWEuWGzWEuit) . (110) 

The above result can be further simplified, if we assume all two-point functions to be 
translation invariant. Then, after a Fourier transformation, the quantities appearing 
in the arguments of the determinants in eq. ( |1U9[ ) can be handled like ordinary 
numbers, which yields: 

SfW = 1 " det (I^hMhV' 2 = l-exp{|lY !„(...)} 



1 + G'iSi 2 G' 2 Si 2 




Gi(k, t)G 2 (k, t)G± 2 (k, i) V 



1 + G^M^M^M) 



where = J d d x in d dimensions. 

Note that for any translation invariant two-point function F here, with 

F(s,y;f) = F(x-y;t) = J ^ e*<«-»> F(k; t) , (112) 

the following constraints, 

F(k;t) = F(-k;t) , and F(k; t) > , (113) 

are implied by the requirements of being symmetric and positive definite. 

Before we continue, we want to draw attention to several remarkable features of 
the above result for the linear entropy, which seem to be of a rather general nature: 

I. Neither the "mean fields" 0i j2 , nor their conjugate momenta 7fi j2 , nor the imaginary 
parts Si j2 of the parton and gluonic environment two-point functions, respectively, 
contribute to the entropy, i.e. to Tr p-p. This result is formally true for arbitrary 
powers, Tr p-p, since the corresponding terms always cancel between successive factors 
of p-p, when evaluating the trace. Thus, we anticipate it to hold for the statistical 
entropy as well (cf. below). 

II. As expected, the entropy vanishes for vanishing correlations between partons and 
gluonic environment, i.e. Gi 2 (t) = Si 2 (t) = (independent subsystems). 

III. If either one of the widths of the parton or environment wave functionals vanishes, 
Gi(t) — > or G 2 {t) — > (cf. eqs. ( |9"5| , |104D , then the entropy vanishes. Thus, if either 
subsystem is constrained to follow essentially the classical equations of motion, then 
there is no entropy production (cf. discussion at the end of Sec. 4.4). 

IV. The above conclusions I. - III. and the result for the linear entropy as given in eq. 
( |1 1 1| ) are model-independent. They hold for any scalar field theory of partons coupled 
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to gluonic environment modes independently of the specific form of the (renormaliz- 
able) interactions. They are, however, obtained here in the TDHF approximation as 
embodied in the wave functional of eq. ( |104| ). 



If we want to learn in detail about the time-evolution of the system, i.e. of the vari- 
ational parameter functions, then we have to specify the interactions. The dynamics 
of the system with the model interactions from Sec. 4.1 will be studied in Sec. 4.4. 

It is worthwhile stressing the importance of parton-environment correlations one 
more time, since these are precisely the quantum correlations discussed after eq. 
(|j~7|) in Sec. 2. Thus, our considerations presently confirm the general idea that 
quantum decoherence in a subsystem, and therewith a measurable entropy production 
therein, is induced by interactions with a dynamically active environment |] The 
parton model serves as a study case par excellence, since confinement dynamically 
separates the total system into an almost directly observable hard parton subsystem 
plus an unobservable non-perturbative environment of soft gluon modes (cf. Sec. 
4.1). 

Finally, we want to present here also an approximate result for the statistical en- 
tropy, which will later be more closely related to observables, the particle multiplicity 
in high-energy reactions, in particular. With the general definition of the statistical 
entropy in eq. ([18]) we find, 



S-p = -^2 c p In c p = (- In c p ) > - In (c p ) 
v 

= -In c p = - ln Tt -pPv = -hi(l-S* n ) , (114) 



where we also used eq. (|19| ) and the inequality results from Jensen's inequality for 
convex functions, which is the basis for many variational principles (see e.g. Ref. |23[|). 
The equality sign holds, if the (sub-) system is in a pure state (vanishing entropy) or if 
all states are equally probable, c p = const, which is the case for a thermalized system 



at asymptotically high temperature (T — > oo), cf. eq. (|23|) for example. Applying 
eq. (|114|) to eq. (|111|) , we obtain: 



o m > l v f d d k ( l-G^G^G^k,^ 

Sv{t) > -~V d J j^y d ^[ 1 + Gl{kjt)G2{kjt)E 2 2M ) ■ (H5) 

The right hand side here should provide a useful approximation for the statistical 
entropy in our scalar model, if we are interested particularly in the initial phase 
of entropy production, i.e. the decoherence process in a strong interaction at high 
energy, starting from a pure quantum state. 

4.4 Effective Action and Equations of Motion 

Our aim in this concluding section is to derive the effective action for the complex 
system of partons and gluonic environment from which we then obtain the relevant 



9 Note that a functional diagonalization of the parton density functional, see eqs. ( |106| ) - (108), 
would amount to an explicit construction of the time-dependent pointer states (see Sec. 2) for the 
parton field. 
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TDHF equations of motion, which describe the time-evolution of the system, the 
partonic subsystem in particular. 

First of all, it helps somewhat to simplify the following calculations, if we rewrite 
the partonic density matrix from eqs. (|107J108|) . We find: 



(0ilM*M> = MMiA (He) 

= *^[0i;t]* Gl [0i;t] e^{Y*[<f>' x ;t}G 2 (t)Y x [(f> x ,;t}} , 

with 

Y x [frt\ = Itf-falGuW-iEnit)] , (117) 

and where the effective gaussian is defined as before, however, with the following 
replacements: 

Nl (t) _ N ± (t) = N x (t)N X2 (t) , 

G?{t) — > G^it) = Gf(t)A(t) , 

Sl ( f ) _> Ei(t) = E 1 (0-|[Si 2 (t)G 2 (t)G' 12 (t)+G' 12 (t)G 2 (t)E 12 (t)] ,(118) 

see eq. fl9~5|) for the definition of the generic and (|110| ) for A(t). (Recall that 
products of two-point functions involve an integration over intermediate coordinates.) 
Not surprisingly, the changes presently induced in the parton density functional by 
the environment bear some distant similarity with what happened to the single-parton 
density matrix in Sec. 3.1, cf. eqs. Q5T| ) - (pT|)- It is also worthwhile to appreciate 
the analogies and differences between the Feynman- Vernon approach for the non- 
relativistic toy model, eqs. (]31f) - (j35D , and the TDHF result for any two coupled 



scalar fields, eqs. ( |116| ) - ( |118| ) above. 

Next, we want to calculate the analogues of eqs. (p8|) - (|100|) to see the changes 
caused by the environment interacting with the scalar parton field: 

(<f>i(x)) = Jv<j>! <h.{x) pv[fa,fa;t} = Mx,t) , (119) 

0(pi{X) J o<p\[x) 

Thus, there is no change in the mean fields here, as expected. However, 

(0i(x)0i(y)) = / £>0l 0i(x)0i(y) p P [0i,0i;t] 



(n x n X2 ) 2 J vfa [Mx) + M^t)}[My) + My,t)} 

■ expj-i^Gr 1 ^-^} 
fa{x,t)fa(y,t)) - 2{N x N X2 f 5 (N x N X2 y 

0Lr 1 [A — IJWfay) 



2 



M^t)Mv,t)} + [A(t) - flWr'GiWk*,*) , (120) 
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where we used eqs. (|97| , |105|) , and eqs. ( |110|) ; we suitably shifted the field to obtain 
the second equation and employed (for M symmetric) 
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5M(x,y) 



det M 



6 



5M(x,y) 



exp{tr In M} = M~\x,y) det M 



(121) 



to obtain the final result. Obviously, the correlations with the environment modify 
the parton two-point function in a non-trivial way. Note that the relevant factor 
[A — B}' 1 could be expanded in a geometric series. We remark that equations which 
are formally identical to (|119|) - (|120[) hold also for the field 2 . Due to the formal 
symmetry of the wave functional, eq. ( |104| ), they are obtained by simply exchanging 
the indices everywhere ("1 <-> 2"), which denote the parton and gluonic background 
fields, respectively. Furthermore, we obtain: 



(idt) 



Vfa V<p 2 #J 2 [0i,0 2 ;t] (id t ) ^i 2 [0i,0 2 ;t] 



rtr 



Sl 2 (t)([01-0l(t)][0 2 -0 2 (t)]) 



(122) 



where we employed eq. ( |120| ). Equation ( |122| ) generalizes eq. ( |100| ) by adding up 
separate contributions from both fields and a new term entirely due to the interaction 
between them. Consequently, a similar structure will arise in the evaluation of the 
total effective action below. We also find: 



h(x) - 4>i(x,t)}[(f)2(x) - (f> 2 (x,t))) 



™12T?A V 

bG 12 {x,y) 
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-l{G 2 (t)G l2 (t)[A(t) - BitT'dit) + [A(t) - B^G^GMG^t)}^ 

= -\A{t) - B(t)]~ 1 G 1 (t)G2(t)G 12 (t)\( x , y ) , (123) 

where the index 5* indicates full symmetrization w.r.t. G\ <-> G 2 and separately 
x <-> y. This symmetry stems from the symmetry of G\ 2 and eq. (|105|) , which 
is frequently used in the course of the calculations here with N£ =det(A — B). 



For translation invariant two-point functions the last equality in (|123|) follows. This 
simplifying assumption will always be made henceforth. 

As a final ingredient for the evaluation of the effective action we need 



i(x) -0i(x,t)] 2 [0 2 (x) -Mx,t)} 2 ) 

( S 2 



(N 1 N 2 Ni 2 ) 2 • 



+ 4 



8Gi 2 (x, x) 2 5G 1 1 (x,x)5G 2 1 (x,x) 



{NtNzNu) 



- 0l][02 - 02] (,)) 2 + ([01 - 0l] 2 ,))([02 - M\ x) ) 



+ 2G\ 2 (x } x) ■ [A — B] G l G 2 Gi 2 \( x ^ x ) + 2GiG 2 Gi 2 \( x ^ x ) ■ [A — B] ~G\G 2 G 



I2\(x,x) 
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+ {G x {x,x) ■ [A - B}- 1 G 2 \ {X , X) + \G x {x,x) ■ [A - B^G^iG?^ 
+\G l GiG^\ (x , x) ■ [A - B]- 2 G x \ {x , x) } + { "1 <-> 2" } , 



(124) 



i.e. an expectation value with all four fields at the same point (the factor 2 in the 
first equation is a symmetry factor like the factor 3 in eq. ( |103j )); we made use 
of eqs. ( |120| , |123| ) to rewrite terms as "disconnected parts" in the first line in the 
second equation and remark that the above result would look even worse without 
a final simplification due to translation invariance (intermediate spatial integrations 
are suppressed as before and the obvious time dependence is omitted). 



To calculate the effective action, eq. (|94]) , for the scalar field theory defined by eqs. 
(^,0) in TDHF approximation we proceed in analogy with Sec. 4.2 and evaluate 
various contributions separately: 



r[tf 12 ] = J dt{^ 12 {t)\[id t -H]\^ 12 {t)) = J dt{{id t ) - (Hi) - (H 2 ) - (V 12 )} . 

(125) 

First of all, (id t ) was evaluated already in eq. (|122p together with eq. (|123|) . Secondly, 
similarly as the corresponding terms (without time-derivatives) in eq. ( |101| ) before, 
we obtain after a lengthy calculation (j = 1, 2): 



d d x <tt 12 (*)|[- 



i 5 2 



+ 1 2 (V<j> j ) 2 +v j (<f> j )]\* 12 (t)) 



(126) 



- tr 



-2Ej[A - B^GjEj - ~Gf x [A + B] + \[A- B]~^G 3 



(2), 



( 4 )/ 



+ Jd?x[ yVy'Wj) ■ [A - BrGjl^ + -fi-'tfj) ■ ([A - BY x G^ x) f } , 
where the same comments as after eq. (|1(J1|) apply and we made use of equations 



analogous to ( |102| , |103|) which hold for the present case. Furthermore, we employed eq. 
( |12(J|) and recalled the definition of Ex i n ec i- (HH), while S 2 follows with the familiar 
exchange "1 <-> 2". We remark that eq. (|126|) correctly yields the corresponding 
terms in eq. ( |101| ) for vanishing correlations (A = 1, B = 0) between the fields </>i 
and 02, i-e. for independent subsystems. 

Next, we have to calculate the contribution (V12) from the mutual interaction 
of the fields. Employing the appropriate Taylor expansion for a general potential, 
which is at most quartic in the fields (cf. also eqs. ( |102| , |103|) ), and keeping only 
non- vanishing terms yields: 



(V12) 



d d x (*i 2 (*M0i,0 2 )|tf 
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Vfa, h) + - 00(02 " U) + ^ (2 ' 2) ((01 " 0l) 2 (02 " 2 ) 2 ) 

+ { ^ (2 ' O) ((0i - 0i) 2 )> + ^ (4 ' O) ((0i - 0i) 4 ) } + { "1 - 2" } , (127) 
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where V^ m ' n ^ = -^j^-^j^V(4>i, 4>2) and we omitted the overall / d d x in the last equa- 
tion. The expectation values arising here have all been calculated before. In particu- 
lar, we recall eq. ( |103| ), which applies for the terms in the last line with ip — > 0i )2 and 
G -1 — > Gi l ' [A — B] after suitably shifting the fields; these terms obviously generate 
additional contributions of a similar (mean field dependent) structure as the poten- 
tial terms oc i>j in eq. ( |126 ). The most interesting terms besides the classical 



contribution in the second to last line in eq. (|127 ) are obtained explicitly by inserting 

eqs. drags. 

This completes our derivation of the Cornwall- Jackiw-Tomboulis type effective ac- 
tion [0] for any system of two coupled scalars fields with at most quartic interactions. 
We implemented Dirac's time-dependent variational principle, eqs. ( p3| , p4|) , with the 



help of a most general gaussian trial wave functional, see eq. (|104|) . Now, combining 
eqs. ( |122| ) - ( |127| ), we obtain the final result (with h's inserted): 



r[*ia] = fdt { £ { TT^ - ±7fJ - i(V^) 2 - Vj fa) 



3=1,2 



+ htr[ Eft - ZLjGft - -Gj [A + B] + -V 2 G 



+ ~ tr [ t l2 G l G 2 G l2 } - Vfa, 2 ) + h^{V^) tr [G 1 G 2 G 12 ] 



V d 4! 



2( tr GiG 2 Gu) 2 + tr G 1 tr G 2 



+2 tr G 12 tr [G x G 2 Gi 2 [A - B}' 2 } + 2 tr [G^G^] tr [G^G^A - B\ 
+{ tr Gi tr G 2 + - tr Gi tr [GiG 2 2 G 12 L4 - B]-' 2 ] 

+ 1 - tr [GiG 2 2 G 1 2 2 ] tr [GM - B]~ 2 ] } + {«!<- 2" } 



128) 



with vf ] = drvjtyj/cty?, Vj n) = rf n 1/(0i,0 2 )/rf0/, and as defined after eq. 

QI271). We also use the abbreviations G i2 = G i2 L4 - B], Gj = [A - B\~ x Gj, and 
V 2 acts on either one of the two formal arguments of Gj] Ex was defined in ( |118| ), 
E 2 follows by "1 <r-> 2", and A,B are given in eqs. (|110| ). We remark that due to 
the assumption of translation invariance of the two-point functions their position in 
a product is irrelevant, if they are connected by integrations over intermediate coor- 
dinates.^ Furthermore, we employed eqs. ( |136| ) from the Appendix and introduced 

10 Implementing translation invariance into the effective action does not interfere with taking one 
further functional derivative as in deriving equations of motion, if proper care is taken of the "zero 
mode" , see the Appendix. 
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the spatial average of a function /, (/) = / d d x f(x)/Vd, with V d = J d d x denoting 
the volume of the region of integration. Of course, all one- and two-point functions 
in eq. ( |128| ) are time- dependent, which we suppressed. 

We conclude this section by deriving the corresponding TDHF equations of motion 
for the ten variational parameter functions entering the effective action in eq. ( |128j ). 
In particular, these equations describe the time-evolution of our parton model field 
theory defined by eqs. (BBlETp together with eqs. 
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SlTj (x, t) 
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X 1 ) 



(1) 
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-'i^f + VP}G 5 {x,x,t) + h^-V^G^G 



r(3) 



d 



(129) 
(130) 



(131) 
(132) 



I2\(x,x\t) , 



where we used eq. (|129|) to indicate that ftj can be easily eliminated from the set of 
equations. Apart from eq. ( pL 3 1| ) and the last term in eq. ( |132[ ), the above equations 
still have a structure similar to the ones which can be obtained from eq. ( |101| ) in the 
case of a single scalar field [[II], [EJ. Note that eq. ( |131| ) can be explicitly written as 
an equation for Gi 2 , eliminating Gj by eq. ( |130| ). 

We remark that the equations following by variation w.r.t. the two-point func- 
tions Ej and Ei 2 , in particular, could be derived by directly applying the respective 
functional derivatives on the effective action, since all contributions arise in this case 
from trace terms in eq. ( |128|) . For the other two-point functions, Gj and Gu, there 
are terms with potential insertions. Here one has to be more careful, taking into 
account that the two-point functions actually depend only on one variable due to 
translation invariance. The relevant formulae are given in the Appendix. 

Then, we obtain another pair of equations which generalize the corresponding one 
for a single scalar field (j' ^ j, j = 1, 2): 



sr 



SGj(x } y } t) 







^j(x,y,t) = [—2T lj T lj + -Gj 2 + -Gj'G^t,^ + Gj'T.j'G^T.u]^^) 



+ 



iv-i«„f> + if»-AI { „<« + vf} tr G J 



(*) 



S(x - y) 
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+ (vW)G r G 



12\(x,y,t) 



A iy(2,2) 
Vw 4 



4 tr [G1G2G12] • G r G 



12 



+ tr Gf ■ [A - B]- 1 + - tr Gj ■ G 2 G 12 F + tr [Gy [A - B}' 1 } ■ G 1 G 2 G 1 



+ tr [G 2 G r G 2 2 ]-G 2 G 2 2 [A-B\ 



tr [GjfA-S] 
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<jyLr 12 



+4 tr Gi 2 • GjG}G 12 + 2 tr [GaG 2 G 12 ] • G,vG 12 [A - £]' 
+2 tr [G1G2G12] • GjfG 12 F + i tr [G.G^G^] ■ [A - £] s 



+ - tr [GjG^G^] • F + tr G,v • (i + [A — £?] 2 ) 



(133) 



where we used formulae from the Appendix and inserted the abbreviations introduced 
after eq. ( |128| ) and F = [1 + GiG2G 1 2 2 ]/[l — GiG2G 1 2 2 ] to write the result as compact 
as possible. For the same reason we collected the arguments from individual terms 
and wrote 1 for a 5-function in the last line. At last, we find: 
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2 tr [G1G2G12] • [i + GiG 2 G 



12J 



+ tr [G 2 G 2 G l2 [A - B\\ ■ G^G^ + tr [GiG 2 G 12 ] • [A - B] 2 

+ tr Gi 2 • G!G 2 (4[/4 - fi]" 1 - 3) + tr [G1G2G12] • (4L4 - F]" 1 - 3) 

+{ tr Gi • G 2 G 12 + tr G x ■ G 2 G 12 (1 + ^F) + tr [GiG 2 2 G 1 2 2 ] • G^ 



+ - tr [Gi[A - B}- 1 } ■ G 2 G 12 [A - B] } + { "1 <- 2" } 



(134) 



where the same comments as after eq. ( |133|) apply; in particular, note that V 2 acts 
on either one argument of the following term, as before. 

( |134| ) we obtained a closed set of ten coupled non-linear 



Thus, with eqs. ( |T2~ 

equations which describe the time-evolution of any theory of two coupled scalar fields 
in terms of one- and two-point Wightman functions (TDHF approximation). They 
are local in time and of first order in time derivatives (see, however, eq. (|132j) for the 
elimination of 7fj and consequently d t (j)j — > d 2 (f>j). 

Generally, one has to expect the usual ultraviolet divergences in terms like tr G(t) = 
J d d x G(x,x) corresponding to expectation values of two fields at coinciding points. 
We assume them to be regulated by a high- momentum cut-off at present. Note that 
by their definition via the Momentum Space Mode Separation in Sec. 4.1 the gluonic 
environment modes are automatically regulated in this respect. Therefore, one could 



36 



essentially take over the separate and only necessary renormalization of the scalar 
parton sector from Refs. |30], |33|] . We plan to come back to this feature in the sequel 
to the present work. There we also have to specify physical initial conditions for the 
set of parton-environment TDHF equations and study the restrictions imposed by 
renormalization. 

In passing we remark that due to translation invariance eqs. ( P-29Q - (|134j) simplify 
considerably after a Fourier transformation. From multiple convolution type spatial 
integrals involving two-point functions there is always only one final momentum inte- 
gration left over, cf. e.g. eqs. ( 1111 ,115) where this was employed before. Incidentally, 



for any translation invariant two-point F function here, with 

F(p;t) = J d d x e~ ip y F(y;t) = J d d x e~ ip - y F(x + \y, x - \y; t) , (135) 

the Fourier transform coincides with the non-covariant Wigner function. Therefore, 
after a straightforward Fourier transformation, our set of equations at the same time 



presents the transport theory |3l] for any two coupled scalar fields in TDHF ap- 
proximation. It is obtained here in the non-covariant form suitable for initial value 
problems and should allow interesting comparisons with other transport equations 
derived in weak-coupling perturbation theory. 

Finally, we point out two special cases of the parton-environment TDHF equa- 
tions, which are of interest in their own right: 

(a) Setting (fij = = 7fj, i.e. no classical mean fields, the set of equations reduces 
to six coupled equations. This may be relevant for a short-lived system with no such 
fields in the initial state which is in some sense dominated by incoherent radiation, 
i.e., if there is no time to develop classical expectation values of the fields. It might 
also apply to a densely (in phase space) populated parton system, if (and only IF) 
the effect of such a "quasi high-temperature" state is to "melt" the classical fields. 

(b) Conversely, if we constrain one part of the complex system to follow the classical 
equations of motion, i.e. to be completely mean field dominated, by setting e.g. G 2 = 
= S 2 , then the set of equations again reduces to only six coupled equations, since the 
dependence on the correlation functions G\2 and £12 is automatically eliminated from 
the effective action, eq. ( |1 28|) . This is in agreement with our conclusions concerning 
entropy production, cf. eqs. ( |111| , |115| ) and point III. in Sec. 4.3, which is identically 
zero here. Obviously, this special case (i.e. the QCD analogue to be worked out) 
cannot be relevant for high-multiplicity events in strong interactions at high energy 
and the parton model as applied to such events, in particular. 

We further discuss these issues in the following section. 



5. Conclusions 

At the outset we apologize to all practitioners and followers of the QCD parton 
model. We did not consider QCD partons and their gluonic environment yet. How- 
ever, our considerations and the basic idea towards a solution of the about 40 years 
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old 11 entropy puzzle" in strong interactions at high energy are motivated by what 
we believe to be important features of QCD. This explains our abuse of the words 
"parton" and "gluonic" environment modes. 

To begin with, throughout this work and in Sec. 2, in particular, we consider the 
von Neumann or statistical entropy defined in terms of the relevant density matrix, 
see e.g. eq. fljjf) , the part of Sec. 3 following point (iii), and especially Sec. 4.3. 
For technical reasons we often employ the linear entropy, eq. flT9|), which we show to 
provide a lower bound on the statistical entropy, see eq. ( |1 14|) . Concerning quantum 



properties of the system to be characterized by the entropy, both definitions are 
equally valid for our purposes. In particular, both measure the impurity of the parton 
density matrix. 

The core of the entropy problem is to try and understand how (an idealized exam- 
ple) two hadronic scattering in-states undergoing a hard interaction, i.e. a quantum 
mechanically pure initial state, can result in a high-multiplicity event corresponding 
to a highly impure (more or less thermal) density matrix on the parton level before 
hadronization. 

Thus, attempts to associate the apparent entropy production with a "coarse grain- 
ing" either in the phase space of the observed system or due to the finite resolution 
in any experimental measurement, which both may be useful to characterize derived 
"macroscopic" aspects of such reactions, seem to miss the point. There one gives up 
from the beginning the possibility to understand on a fundamental dynamical level 
how a complex pure-state quantum system can produce classical behaviour, i.e. an 
impure density (sub-)matrix with decoherence of associated (parton) wave functions 
and entropy production. 

In Sec. 2 we introduce a convenient general framework in terms of Schmidt and 
pointer states to analyze the consequences of a dynamical separation of a complex 
(possibly strongly interacting) system into an observable ( "partonic" ) subsystem plus 
unobservable ("gluonic") environment modes [|7|, §, ||, pT| . Entropy production due 
to environment-induced quantum decoherence in the observable subsystem arises nat- 
urally, IF such a separation of the total "closed" system into an "open" subsystem 
and its environment is dynamically realized. There is no guarantee for this to happen 
and the existence of the associated (almost) classically behaving (partonic) pointer 
states, i.e. states the quantum superpositions of which are dynamically suppressed, 
has to be viewed as a particular feature of the system. 

Such a major miracle seems to be installed in QCD [p9fl : for example, the fact 



that deep-inelastic scattering can be described accurately in terms of a hard scattering 
cross section and structure functions. The latter are decoherent one-particle probability 
distributions. Therefore, some "secret agent" has to effect this extremely efficient 
decoherence process, once a parton is knocked out of its coherent initial state wave 
function. Our point of view here does not conflict with the standard parton picture, 
in particular, the applicability of plane-wave states representing partons entering or 
leaving Feynman diagrams.^ Any basis should do. However, the miracle consists 
in the fact that typically the initial state of a hadronic scattering reaction can be 
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I thank R. Baier for insisting with his questions about this. 
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described by structure functions which correspond to diagonal density matrices, i.e. 
without interference terms (which always can be achieved formally by the Schmidt 
decomposition procedure, however, only at one instant of time, cf. Sec. 2). 

In terms of the split hamiltonian in eq. (plj), the bound partonic initial state 
wave function diagonalizes H-p + H-pg. Once it is perturbed by an external (e.g. 
electromagnetic) interaction, the interaction with the gluonic environment starts to 
"rattle and shake" the perturbed wave function so strongly that it decays very fast 
into a decoherent superposition of states as reflected in a structure function. 

We conjecture that gluonic modes (neglecting quarks for now) which are close 
to constant field configurations form an essentially unobservable sector of the total 
Hilbert space of a QCD system. The associated Momentum Space Mode Separation 
is discussed in more detail in Sec. 4.1. They act as an active environment on the 
observable^ partonic sector, to which they are strongly coupled via low-|Q| 2 processes. 
Presently we explore this picture in two simple models. 

In Sec. 3 we study a non-relativistic parton coupled to a gluonic environment in 
close analogy to an electron coupled to the quantized electromagnetic field [^7]. |2"Bf . 
However, we deliberately change the spectral density of the environmental oscilla- 
tors to be dominated in the infrared. We calculate the time-dependent parton den- 
sity matrix by integrating out the environmental degrees of freedom exactly with 
the Feynman- Vernon influence functional technique |23|, |26 ]. In a short-time strong- 



coupling limit, which has not been of interest in studies of quantum Brownian motion 
so far f2"6fl , we find the following results analytically (cf. fig. 1 and points (i) - (iii) 



in Sec. 3.2): 

Gaussian parton wave packets experience friction and localization and their coherent 
superpositions decohere. 

All effects are governed by a non-perturbative time-scale f^ 1 <C lfm/c, eq. (0). 
Thus, a parton following the classical trajectory with the center of its wave packet is 
slowly decelerated, which can be interpreted as being due to the emission of infrared 
gluons, i.e. the excitation of environment modes in this model. More surprising is 
the localization effect, i.e., the quantum mechanical spreading of the wave packet can 
be suppressed or even reversed depending on the choice of parameters. The initial 
wave function obtains an almost soliton like character through the strong interaction 
with the environment. Finally, considering the superposition of wave packets, we 
find that their decoherence actually happens in the short-time strong-coupling limit. 
Consequently, in accordance with the general considerations of Sec. 2, we "see" and 
can calculate how entropy is produced in this toy model on a short time scale. 

All the above effects, if recovered in QCD eventually, seem to be very important 
e.g. for a justification of the classical cascade approach to multiple parton scattering 
|18|| . Then, a parton propagating in space-time behaves essentially like a classical 



particle in between perturbative hard scatterings. 

In Sec. 4 we set up a fairly general model of two coupled scalar fields representing 
partons and gluonic environment modes, respectively (see Sec. 4.1 for details). In 
this case, due to the non-linear interactions, the environment degrees of freedom 
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In the sense of parton-hadron duality or deep-inelastic scattering. 
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cannot be integrated out as before. Therefore, our main concern here is to formulate 
a tractable approach to the time-evolution problem in this quantum field theory, 
which is not restricted to perturbatively small couplings, in particular. We employ 



Dirac's variational principle [|T^, 16|, see Sec. 4.2, and describe the complex system 
in terms of its ten one- and two-point Wightman functions (TDHF approximation). 
This corresponds to the most general gaussian wave functional, eq. ( |104| ), in the 
field theory Schrodinger picture. In Sec. 4.4 we derive the related Cornwall- Jackiw- 



Tomboulis type effective action [[34], eq. (|128|) , governing the dynamics and obtain 
the equations of motion from it. 

An attractive feature of having an explicit wave functional in hand is that it allows 
to calculate the partonic density functional here quite generally, eqs. (|106|) - (|108|) 
or eqs. ( |116| ) - ( |118| ), without having to solve the complicated dynamics first of all. 
Thus, we obtain in Sec. 4.3 a model-independent result for the entropy in terms of 



two-point functions, eqs. ( |lll|jll5| ), which causes several interesting observations, see 



points I. - IV. there, which are presumably of a general nature: 

1. If there are no quantum correlations between the partons and the gluonic environ- 
ment, then there is no parton decoherence and no entropy production. Again, this is 
in accordance with the general considerations in Sec. 2. 

2. If one assumes a strictly classical gluonic environment, i.e. following the classical 
equations of motion, then there is no parton decoherence and no entropy production. 

This latter observation is in interesting contrast to recent attempts to explain 
entropy production and thermalization in QCD as being due to the non-perturbative 
chaotic dynamics of purely classical Yang-Mills fields, see Ref. [^] and references 
therein. There the above mentioned "coarse graining" of the classical phase space 
covered by an ensemble of identical systems^ is necessary to deduce a classical en- 
tropy. The results are very suggestive as compared to high-temperature field theory. 
However, in both cases decoherence is put in by hand. Their relevance for the von 
Neumann entropy considered here and the most important dynamical decoherence 
phenomena in a quantum parton system, in particular, seems to be unclear at present. 

Our results in Sec. 4 will be the starting point of the sequel to this work, where 
we plan to study physical initial conditions for the parton-environment system and 
their relation to model structure functions. Then, of course, the time-evolution of the 
system following the equations of motion, eqs. (|129| ) - ( 134 ), will be most interest- 



ing to consider. Any solutions will provide the explicit time-dependence of entropy 
production according to eqs. ( |111| , |115| ), especially for early stages of a reaction. 

Let us summarize our present point of view by saying that a parton appears like 
a parton only because it feels the gluonic environment which manifests its strong 
non-perturbative interactions caused by the running coupling on a short time scale 
(^C lfm/c) through the induced decoherence properties of partons and the associated 
observable entropy production. 

Finally, it seems worthwhile to mention a few other potential applications of our 

13 Of course, this does not represent a pure quantum state and, thus, the resulting "thermalization" 
seems to be a very plausible consequence of the mixing property of a strongly chaotic system, i.e. 
the filamentation of the phase space occupied by the ensemble. 
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results from Sec. 4. They provide a first step to study the analogue of quantum 
Brownian motion []26|| in the context of field theory. In particular, one may study 
the time-dependence of a Higgs model type phase transition under the influence of 
the interaction with a perturbing environment field. The possibility of such a non- 
equilibrium phase transition is contained already in our simplest specification of the 
interactions in eqs. (S8|) . It has been another long-standing problem to study changes 
of signals, which are thought to characterize an equilibrium phase transition, when 
they come from an out-of-equilibrium system such as the early stage of a high-energy 
nuclear collision. Recently, following the suggestion of the possible formation of a dis- 
ordered chiral condensate in relativistic heavy-ion collisions, it has been realized that 



the condensate field actually presents an open system, see Refs. |36[ and numerous 
references therein, to which our methods can be applied. 
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Appendix 

We list several formulae which are useful for the calculations in Sec. 4.4. 14 They 
are valid for any translation invariant symmetric two-point functions X, Y, Z, and 
any integrable function V (in any dimension): 

I = J dxV(x)-XYZ\ (x>x) 

= J dx dx\dx2 V{x)X[x — x\)Y[x\ — X2)Z[x2 — x) 

= J dx V(x) ■ XZY\(x tX ) [ or any other permutation of XYZ ] 

dxV(x) tr [XYZ]/vol , (136) 



where vol = f dy denotes the volume of the region of integration. Then, functional 
derivatives w.r.t. one of the two-point functions can be calculated in two ways em- 
ploying the trace formula in (|136|) : 



5 , r 6 



[ dx' V(x') —— tr \XYZ]lvol 
J oZ(z) 



5Z{z) J v ' 5Z{z 



14 C. Wetterich helped with a remark here at the right moment. 
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SZ(x,y) 

Dividing out the constant factor, one obtains: 



J dx V(x') XY\{x tV ) [for z = x — y ] 

I dx ' v ^jzL^ aixYz] - (m) 



tr [XYZ] = vol ■ — tr [XYZ] , (138) 



5Z(x-y) ' SZ(x,y) 

which clearly exhibits the "zero mode" factor connecting the functional derivative 
which takes translational invariance into account with the one which does not. 



Finally, with the definitions of A, B and Gj, G\i given in eqs. (|110|) and after eq. 
|128Q , respectively, one finds the following simple results: 



(139) 

where the permutability of translation invariant symmetric two-point functions in a 
product, e.g. XYZ\^ y \ = XZY\/ X ^ = ... (cf. also eq. ( |136| )), is crucial. Similar 
useful relations are easily obtained for any power of [A — B] multiplying Gj or G12. 



Figure captions: 

Fig. 1: Deceleration and localization of a non-relativistic parton in a gluonic envi- 
ronment. Shown are the unique deceleration v/v , eq. fl66|) , and the relative change 
of the width (w — /?)//?, eq. (|69|), of the wave packet as a function of t + = f + t. We 
set a< = 0.1a; see eqs. ( ffO|) for the definitions of the parameters. For comparison 
free particle results (dashed curves) are given. 

Fig. 2: Qualitative picture of Momentum Space Mode Separation: "Unobservable" 
gluonic environment modes are confined to the inner region II (roughly fc 2 < and 
k 2 < Aq CD ). Partons live in the outer region I. 
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